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Abstract:  This  report,  the  first  of  a  series,  describes  the  coastal  modeling 
package  TWAVE  (Typhoon  and  WAVE)  developed  under  the  Surge  and 
Wave  Island  Modeling  Studies  (SWIMS)  project  capable  of  modeling 
tropical  cyclone  winds,  waves,  storm  surge,  and  coastal  inundation.  The 
TWAVE  flood  and  inundation  package  is  a  practice-oriented  engineering 
modeling  tool.  The  first  version  of  TWAVE  is  described  in  this  report  and 
will  be  updated  based  on  user  feedback  in  subsequent  reports  and  as  new 
modeling  capabilities  become  available.  The  objective  of  this  research  is  to 
develop  a  system  with  multiple  levels  of  complexity  that  is  suitable  for 
applications  such  as  hindcast,  forecast,  and  hypothetical  storms.  TWAVE 
is  a  personal  computer  based  modeling  system  and  uses  Microsoft  Excel® 
to  organize  and  visualize  input  and  output  data.  The  package  includes 
three  wind  models  coupled  to  a  deepwater  spectral  wave  model  and  a 
parametric  wave  model  with  several  options  for  calculating  reef  top  wave 
heights,  wave  setup,  and  runup.  These  calculations  range  from  simple 
empirical  formulations  to  a  fully  nonlinear  Boussinesq  wave  model. 

TWAVE  modeling  system  is  currently  configured  for  the  U.S.  Territory  of 
Guam  and  the  island  of  Kauai,  HI,  but  can  easily  be  configured  for  other 
locations.  A  full  description  of  the  TWAVE  package  of  programs  and  input, 
and  output  are  described  in  this  report  to  assist  users  who  interested  in 
adopting  the  system  to  their  specific  needs.  TWAVE  is  a  modular  system, 
and  this  makes  it  possible  to  add  new  components  and  replace  existing 
modules  as  necessary.  TWAVE  modeling  package  is  validated  using 
measured  offshore  winds  and  waves,  and  coastal  inundation  records  for 
Hurricane  Iniki  at  Poipu  beach  on  the  island  of  Kauai,  HI.  Modeling 
estimates  generally  compare  reasonably  well  with  measurements.  The 
calculated  offshore  winds,  waves,  and  coastal  inundation  estimates 
compare  favorably  to  data.  Although  there  were  no  data  for  validation  on 
the  island  of  Guam,  TWAVE  capabilities  are  also  demonstrated  for 
Typhoon  Russ. 


DISCLAIMER:  The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication,  or  promotional  purposes. 
Citation  of  trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the  use  of  such  commercial 
products.  All  product  names  and  trademarks  cited  are  the  property  of  their  respective  owners.  The  findings  of  this 
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Preface 

The  Surge  and  Wave  Island  Modeling  Studies  (SWIMS)  of  the  Coastal 
Field  Data  Collection  Program  is  developing  a  coastal  modeling  package 
called  TWAVE  for  tropical  environments  such  as  the  Hawaiian  Islands  and 
territory  of  Guam.  TWAVE  consists  of  a  set  of  numerical  models  for  island 
coastal  inundation,  and  is  intended  as  a  predictive  tool  for  the  U.S.  Federal 
and  State  civil  defense  agencies  and  the  U.S.  Army  Corps  of  Engineers 
District  Offices.  This  report  describes  the  theory  and  implementation  of 
TWAVE  and  contains  examples  that  demonstrate  its  purposes.  The  model¬ 
ing  package  components  are  described,  step-by-step  instructions  are  pro¬ 
vided  on  how  to  run  each  numerical  model  component  of  TWAVE  and 
view  modeling  results. 

This  report  was  prepared  by  Dr.  Zeki  Demirbilek  and  Alejandro  Sanchez, 
Harbors,  Entrances,  and  Structures  Branch,  and  Dr.  Jane  M.  Smith, 
Coastal  Processes  Branch,  Coastal  and  Hydraulics  Laboratory  (CHL), 

U.S.  Army  Engineer  Research  and  Development  Center  (ERDC), 

Vicksburg,  MS.  J.  Holley  Messing,  Coastal  Engineering  Branch,  Navigation 
Division,  CHL,  typed  the  equations  and  format-edited  the  draft  report. 
Work  at  CHL  was  performed  under  the  general  supervision  of  Jose 
Sanchez  (former  Chief),  and  Jackie  Pettway  (current  Chief),  Harbors, 
Entrances,  and  Structures  Branch;  Ty  Wamsley,  Chief,  Coastal  Processes 
Branch;  Dr.  Rose  Kress,  Chief,  Navigation  Division;  and  Bruce  Ebersole, 
Chief,  Flood  and  Storm  Protection  Division.  Dr.  William  D.  Martin  was 
Deputy  Director,  CHL,  and  Thomas  W.  Richardson  was  Director,  CHL. 

COL  Gary  E.  Johnston  was  Commander  and  Executive  Director  of  ERDC. 
Dr.  James  R.  Houston  was  Director. 
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1  Introduction 

Overview 

Tropical  islands  may  encounter  strong  wind  events  due  to  tropical 
cyclones  which  cause  severe  damage  and  life-threatening  conditions.  The 
capability  to  quickly  forecast  the  coastal  wind  and  wave  conditions  for  esti¬ 
mating  the  resulting  coastal  inundation  is  crucial  for  planners  and 
decision-makers  responsible  for  early  warnings,  evacuations,  and  relief 
procedures.  Protective  coastal  measures  are  commonly  based  on  historical 
records  of  the  measured  surge  levels  and  waves  and/or  numerical  esti¬ 
mates  of  hypothetical  cyclones  representative  of  historical  events.  One  of 
the  challenges  in  estimating  storm-induced  coastal  inundation  is  the  reli¬ 
ability  of  predicting  the  nearshore  waves.  This  requires  using  nonlinear 
wave  models  such  as  Boussinesq-type  fully  nonlinear  models  to  obtain 
reliable  estimate  of  waves,  wave  setup  and  wave  runup  on  coastal  areas  of 
islands.  Because  these  two-dimensional  (2D)  models  are  computationally 
demanding,  an  alternative  is  to  use  simpler  one-dimensional  (lD)  models 
to  simulate  wave  transformation  over  several  transects  along  the  coast¬ 
lines  that  are  of  concern  for  flood  inundation.  Nwogu  (2006)  obtained 
consistent  wave  runup  using  lD  and  2D  Boussinesq  wave  model  simula¬ 
tions  of  nearshore  waves  during  Hurricane  Iniki  along  the  coast  in  Kauai, 
HI.  The  lD  models  have  been  used  extensively  and  successfully  in  calcu¬ 
lating  waves  in  the  surf  zone  (e.g.,  Gerritsen  1980;  Thornton  and  Guza 
1983;  Stive  and  De  Vriend  1994;  Massel  and  Gourlay  2000;  Massel  and 
Brinkman  2001;  Gourlay  2005). 

Over  the  past  decade,  the  Coastal  and  Hydraulics  Laboratory  (CHL)  of  the 
U.S.  Army  Engineer  Research  and  Development  Center  (ERDC),  has 
developed  and  improved  numerical  modeling  procedures  for  determining 
coastal  waves  and  inundation  for  tropical  storms  (e.g.,  Luettich  et  al.  1992; 
Smith  1993;  Thompson  and  Cardone  1996;  Smith  et  al.  2001;  Thompson 
and  Scheffner  2002;  Thompson  2005;  Demirbilek  et  al.  2007;  Demirbilek 
and  Nwogu  2007;  Demirbilek  et  al.  2008).  The  Surge  and  Wave  Island 
Modeling  Studies  (SWIMS)  project  of  the  U.S.  Army  Corps  of  Engineers 
(USACE)  is  developing  predictive  capabilities  for  island  flood  estimates. 
One  of  these  tools  is  a  simple,  first-order  integrated  wave  modeling  pack¬ 
age,  called  TWAVE  (Sanchez  et  al.  2007).  TWAVE  is  intended  for  the 
U.S.  civil  defense  agencies  and  USACE  District  Offices,  providing  users  a 
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practice-oriented  package  which  is  user-friendly,  flexible,  robust,  and  suit¬ 
able  for  the  feasibility  and  reconnaissance  type  of  engineering  planning 
studies.  This  is  a  manual  for  the  TWAVE  model,  not  a  full  assessment  of 
the  model’s  abilities,  which  is  a  task  that  is  just  beginning,  now  that  the 
model  is  operational.  The  TWAVE  modeling  system  is  expected  to  evolve 
and  be  improved  as  it  is  used  in  engineering  projects.  TWAVE  will  be 
calibrated  with  additional  field  observations  in  progress.  Other  numerical 
models  and  different  engineering  methods  may  be  incorporated  into 
TWAVE  as  these  validations  warrant. 

The  missions  and  areas  of  responsibility  for  the  USACE  include  many 
islands  in  the  tropical  Pacific  Ocean,  Atlantic  Ocean,  and  Caribbean  Sea. 
These  islands  typically  have  narrow  coasts,  rugged  interiors,  and  coastal 
roads  and  communities  that  are  a  vital  part  of  the  islands  economy,  trans¬ 
portation,  and  safety.  Many  island  coasts  are  exposed  to  long-period  swells 
from  extratropical  storms  generated  in  both  the  northern  and  Southern 
Hemisphere,  and  are  occasionally  impacted  by  powerful  tropical  storms. 
Despite  the  presence  of  protective  coral  reefs,  elevated  water  levels  caused 
by  wave  setup,  winds,  and  low  atmospheric  pressure  during  storms  cause 
waves  to  run-up  on  shore  and  cause  damage  to  island  roads  and 
communities. 

Previous  studies  employed  a  stream  of  up  to  ten  different  engineering 
models.  These  Fortran  programs  had  to  be  run  sequentially,  and  results 
were  visualized  using  external  commercial  plotting  software  packages.  The 
stream  of  programs  begins  with  calculating  wind  fields  and  ends  with  cal¬ 
culation  of  the  water  levels  and  maximum  runup  levels  at  user-selected 
transects  along  a  section  of  the  islands  coast.  The  set  of  stand-alone 
models  has  been  integrated  into  a  user-friendly  island  flood  and  inunda¬ 
tion  modeling  package  called  “TWAVE,”  an  acronym  for  Typhoon  and 
WAVE  modeling  package.  The  TWAVE  system  is  implemented  in  a 
Microsoft  Excel®  workbook  that  serves  as  the  graphical-user-interface 
(GUI)  for  data  management,  program  execution,  and  visualization  and 
analyses  of  results.  Microsoft  Excel  was  chosen  because  it  provides  a  con¬ 
venient  interface  for  the  package,  readily  available  on  desk-top  computers, 
and  is  commonly  used  by  potential  users  of  the  TWAVE  from  Federal, 
State,  and  local  agencies.  This  gives  users  the  flexibility  to  edit  existing 
plots  and  adapt  graphics  to  fit  their  special  needs. 
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The  TWAVE  flood  and  inundation  package  is  a  practice-oriented  engineer¬ 
ing  modeling  tool.  The  first  version  of  TWAVE  described  in  this  report  is 
expected  to  evolve  based  on  its  utility  for  Federal,  State,  and  local  agencies. 
The  objective  is  to  provide  planners,  decision-makers,  and  engineers  with 
a  predictive  modeling  tool  by  making  it  available  in  the  simplest  form  pos¬ 
sible.  TWAVE  provides  users  the  power  of  present  numerical  modeling 
technology,  places  it  directly  in  their  hands,  and  allows  them  to  apply  it  to 
real-world  scenarios.  This  framework  allows  model  developers  to  refine 
and  expand  the  package  into  a  more  comprehensive  tool  for  future  needs. 
Although  the  first  release  of  TWAVE  has  been  configured  for  the  U.S.  Ter¬ 
ritory  of  Guam  and  the  island  of  Kauai,  HI,  the  modeling  package  is  gen¬ 
eral  and  can  be  adapted  for  application  to  other  islands. 

The  water  depth,  wave  sheltering,  and  shoreline  orientation  are  needed  for 
calculation  of  the  nearshore  waves  in  TWAVE.  Wave  runup  and  inunda¬ 
tion  calculations  require  cross-shore  bathymetric  profiles.  Users  can 
define  an  arbitrary  number  of  nearshore  stations  and  runup  locations  by 
providing  these  required  data.  Future  versions  of  TWAVE  will  include  a 
database  of  pre-computed  cross-shore  profiles  and  enhanced  numerical 
modeling  features  and  capabilities  as  options  for  these  calculations. 
Although  TWAVE  is  intended  for  modeling  storm  waves  for  flooding  and 
island  inundation  works,  it  may  also  be  used  for  the  transformation  of 
swell  and  extratropical  storm  waves  in  the  nearshore  region  for  other 
engineering  projects. 

Modeling  package 

The  modeling  approach  in  TWAVE  is  a  multi-level  approach,  in  which 
users  have  the  option  of  running  different  models  with  varying  accuracy, 
resolution,  and  computational  time.  Figure  l  shows  a  flowchart  of  the 
TWAVE  modeling  system  and  components  as  described  in  this  report 
(blue  text).  The  components  that  will  be  added  in  future  releases  are 
shown  in  red.  Many  of  the  modeling  components  in  TWAVE  are  similar  to 
those  described  by  Thompson  and  Scheffner  (2002)  and  Militello  et  al. 
(2003).  However,  in  this  first  release  of  TWAVE,  the  storm  surge  is 
approximated  by  a  barometric  tide  calculated  based  on  the  atmospheric 
pressure.  Such  a  simple  tide  calculation  method  eliminates  the  need  of 
using  a  computationally-intensive  circulation  model.  This  method  is 
considered  adequate  for  most  tropical  islands  that  lack  wide  continental 
shelves  that  can  alter  the  storm  surge  signal.  Subsequent  versions  of 
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Figure  1.  TWAVE  modeling  approach. 


TWAVE  will  include  circulation  models  for  calculating  the  combined 
storm  surge  and  astronomical  tides.  Parametric  wind  and  wave  models  are 
used  for  fast  computations  as  required  within  the  TWAVE  system. 

More  emphasis  is  given  to  detailed  modeling  of  nearshore  wave  processes 
that  determine  the  extent  of  island  flooding  and  inundation.  Consequently, 
a  wave  energy  flux  model  and  an  advanced  lD  Boussinesq  model  are 
included  in  the  TWAVE  system  for  reliable  estimate  of  nearshore  wave 
transformation,  setup,  and  runup.  TWAVE  has  input/output  (I/O)  files 
associated  with  different  numerical  models  used  in  the  package.  Users  do 
not  need  to  be  concerned  about  most  of  these  files,  but  the  names  and 
their  connectivity  with  other  programs  are  described  in  this  report.  This 
information  should  help  users  to  independently  run  specific  components 
of  TWAVE  such  as  the  tide  model  only  if  so  desired.  An  overview  of  the 
TWAVE  main  modeling  components  is  provided  in  this  report,  including  a 
description  of  the  individual  programs,  processing  subroutines,  and  input 
and  output.  All  input  parameters  and  data  used  in  TWAVE  are  in  SI  units, 
and  times  are  the  Greenwich  Mean  Time  (GMT). 
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The  layout  of  this  report  is  as  follows.  Wind  modeling  is  described  in 
Chapter  2.  The  astronomical  tide  database  used  in  TWAVE  is  reviewed  in 
Chapter  3.  Wave  models  are  presented  in  Chapter  4.  The  TWAVE  package 
including  its  Excel®  interface  and  I/O  files  are  described  in  Chapter  5.  The 
procedure  for  running  TWAVE  is  described  in  Chapter  6.  Two  example 
applications  are  provided  in  Chapter  7.  The  summary  and  conclusions  are 
discussed  in  Chapter  8.  References  section  provides  additional  informa¬ 
tion  as  referenced  in  the  body  of  this  report. 
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2  Regional  Wind  Models 

Wind  and  atmospheric  pressure 

Three  methods  are  available  in  TWAVE  for  calculating  the  time-history  of 
wind  and  atmospheric  pressure.  The  first  method  uses  the  Planetary 
Boundary  Layer  (PBL)  model  and  the  other  two  use  the  parametric 
Holland  and  Rankin  Vortex  models.  The  key  features  of  each  model  are 
briefly  discussed  next. 

In  the  PBL  model,  the  effect  of  marine  boundary  layer  is  considered  by 
relating  the  surface  stress  to  the  geostrophic  flow  over  the  ocean  surface. 
The  PBL  as  used  in  TWAVE,  takes  a  synoptic-scale  pressure  field  to  calcu¬ 
late  the  surface  geostrophic  flow  from  it,  and  then  adjusts  (corrects)  it  for 
the  curvature  effects  of  isobars  (gradient  wind  effects).  To  account  for 
stratification  above  the  sea  surface,  the  air-sea  temperature  differences, 
surface  roughness,  and  moisture  due  to  wave  generation  by  wind  stress 
have  been  included  in  the  PBL  governing  equations  (Cardone  et  al.  1992; 
Thompson  and  Cardone  1996).  A  brief  summary  of  the  PBL  model  follows. 

When  less  accuracy  is  required  (e.g.,  simulating  hypothetical  storms  in 
engineering  applications  that  do  not  require  high  accuracy),  parametric 
hurricane  models  offer  a  simpler  and  faster  approach  to  estimate  winds 
and  pressure  fields  generated  by  tropical  storms.  The  simulation  time  for 
the  PBL  model  is  several  minutes,  and  additionally  PBL  requires  pre-  and 
post-processing,  while  the  run  time  for  the  parametric  models  is  seconds 
and  require  no  pre-  or  post-processing.  The  two  parametric  wind  models 
used  in  TWAVE  have  been  shown  to  provide  reasonably  accurate  wind  and 
pressure  estimates  (Phadke  et  al.  2003). 

PBL  model 

The  PBL  model  used  in  TWAVE  is  based  on  the  vortex  model  of  Chow 
(1971)  as  modified  by  Cardone  et  al.  (1992)  and  later  Thompson  and 
Cardone  (1996).  The  model  simulates  wind  based  on  the  vertically 
averaged  primitive  equations  of  horizontal  momentum  in  a  moving 
coordinate  system  as 
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dV 

dt 


fKx[v  -Vg 


=Vp  +  V.[^„VK]-% 


v+v„ 


v+v„ 


(1) 


where: 

V  =  horizontal  wind  velocity  vector  relative  to  the  storm  center 
K  =  unit  vector 

^9  =  constant  geostrophic  velocity  vector  relative  to  the  storm 
center 

p  =  mean  air  density 
p  =  atmospheric  pressure 
Kh  =  horizontal  eddy  viscosity 
Cd  =  drag  coefficient 
h  =  elevation  of  the  PBL  model 
Vc  =  velocity  vector  of  the  storm. 


The  vertical  advection  of  momentum  and  the  shear  stress  at  the  top  of  the 
PBL  are  assumed  to  be  small  and  ignored.  The  model  includes  parameter- 
izations  for  momentum,  heat  and  moisture  fluxes,  and  a  formulation  for 
surface  roughness  and  drag.  The  pressure  field  is  calculated  using  the 
exponential  pressure  law, 

p  =  pc+APexp(-Rmax/r)  (2) 


where  r  is  the  radial  distance  and  Rmax  is  the  radius  of  maximum  winds. 
The  PBL  model  implements  five  nested  sub-grids  and  a  moving  coordinate 
system  so  that  the  center  of  the  storm  is  always  in  the  center  of  the 
domain.  The  model  has  been  shown  to  accurately  model  surface  wind 
speeds  and  directions  for  tropical  storms  over  open  water  (Thompson  and 
Cardone  1996).  For  Guam  and  Hawaii,  the  largest  grid  is  configured  to 
cover  a  region  of  8  by  8  deg  (890  by  890  km)  with  a  5-min  (9.3  km)  reso¬ 
lution.  The  five  sub-grids  have  12  nodes  in  the  x  and  y  directions  with 
varying  resolutions. 

Holland  model 

This  parametric  model  provides  estimates  of  the  wind  and  pressure  based 
on  the  maximum  sustained  wind  speed,  the  radius  of  maximum  wind 
speed,  and  a  peakedness  parameter.  The  advantage  of  this  model  is  that  it 
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is  fast  and  gives  reasonable  estimates  of  both  wind  and  pressure.  Wind 
velocities  are  calculated  from  the  pressure  distribution  using  a  pressure 
gradient  wind  balance  (Holland  1980)  given  by 


V 


r_j_ 

2 


(3) 


where/is  the  Coriolis  parameter  defined  as/=  2Q  sin  9,  where  Q  is  the 
earth’s  angular  rotation  in  radian/sec  and  9  is  the  latitude  in  degrees.  Vc  is 
the  cyclo-strophic  wind  speed  in  m/sec  obtained  from  the  pressure  field  as 


V  = 


R  B— exp (-R  B 

P 


1/2 


(4) 


where: 


B 

R 


R 


max 


p 

p 

A P  =  pa~  Pc 


peakedness  scaling  parameter 
normalized  radial  distance  (JR  =  r/Rmax) 
radius  of  maximum  winds 
pressure  at  radius  r 
mean  air  density 

pressure  deficit  where  pa  is  the  ambient  pressure  and  pc 
the  central  pressure. 


The  radial  variation  of  atmospheric  pressure  is  given  by 

p  =  pc  +  APexp(-R~B)  (5) 


Note  that  in  Equation  5  the  wind  speed  at  R  =  1  only  depends  on  B  and  the 
pressure  deficit.  Therefore,  if  both  the  maximum  wind  speed  and  central 
pressure  are  available,  a  reasonable  estimate  of  the  scaling  parameter  can 
be  obtained  from  Equation  5  by  setting  R  =  1  (Xie  et  al.  2006)  as 


B  ^Uj 

A P 


(6) 


When  only  the  central  pressure  is  available,  an  estimate  of  the  scaling 
parameter  can  be  obtained  from  the  empirical  relation  of  Harper  and 
Holland  (1999)  as  an  alternate  to  Equation  6.  The  estimate  is  given  by 
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d  —900  „ 

B  =  2-^ - for  1<B  <  2.5 

160 


(7) 


See  the  Control  File  discussion  in  Chapter  5  for  additional  information  for 
specifying  a  constant  scaling  parameter  or  variable  to  be  used  in 
Equations  6  or  7. 

Modified  Rankine  vortex 

In  the  modified  parametric  Rankine  vortex  model  (Rankine  1901),  the  dis¬ 
tribution  of  the  wind  speed  is  related  to  the  normalized  radial  distance  R 
and  the  shape  factor  X  as  follows: 


U^R*  for  K  <  1 
U^R-X  for  if  >  1 


(8) 


where  UmSLX  is  the  maximum  surface  wind  speed  which  occurs  by  definition 
at  R  =  1  and  R  =  r/Rmax.  Hughes  (1952)  has  shown  that  the  empirical  shape 
parameter  X  (exponent  of  R )  ranges  from  0.4  to  0.6.  The  shape  parameter 
in  Equation  8  is  different  from  that  of  the  Holland  model  (Equation  4).  In 
TWAVE,  the  default  value  of  X  =  0.5  is  used,  and  this  value  can  be  changed 
by  the  user  in  the  Control  File  using  the  variable  xRankine.  This  model  has 
been  previously  shown  to  provide  reasonable  wind  estimates  for  the 
Hawaiian  Islands  (Phadke  et  al.  2003). 

Wind  corrections  for  forward  motion  and  inflow 

In  the  Northern  Hemisphere,  tropical  storms  have  stronger  winds  on  the 
right  side  and  weaker  winds  on  the  left  side  of  the  forward  direction  of  the 
storm.  In  order  to  simulate  this  asymmetric  wind  pattern,  a  correction  was 
added  to  the  wind  speed  (Georgiou  1985)  which  takes  into  account  the 
forward  velocity  of  the  storm  as 

VF  =  0.51^  sin  (3  (9) 

where  Vf  is  the  hurricane  forward  velocity  and  (3  is  the  direction  of  any 
location  (each  isobar’s  tangent  angle)  with  respect  to  the  hurricane  trans¬ 
lation  direction. 
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In  the  parametric  wind  models,  the  wind  direction  in  the  Northern  Hem¬ 
isphere  has  a  counter-clockwise  rotation  and  is  approximated  by  a  tangent 
to  the  isobars.  This  requires  the  wind  direction  to  be  corrected  for  the  wind 
isobars  that  are  close  to  the  center  of  the  storm  (core  flow  or  inner  flow 
region).  An  appropriate  correction  was  suggested  by  Bretschneider  (1972) 
by  subtracting  the  following  correction  from  each  isobar’s  tangent  angle: 


(3  = 


10° 


1  + 


R 

20°+25° 

25°, 


max  ) 

V 


R„ 


0  <  r  <  R„ 


,Rmax<r<1.2Rn 


R  >r 

max  — 


(10) 


Parameterizations  of  storm  characteristics 

The  simulation  of  wind  and  pressure  fields  in  TWAVE  requires  specifying 
storm  track  coordinates  in  degrees  north  and  east,  and  l-min  maximum 
sustained  wind  speed  in  knots.  The  storm  central  pressure  (mbar)  and 
radius  of  maximum  wind  speed  (km)  are  optional.  If  the  central  pressure 
is  not  available  (denoted  by  -99),  the  central  pressure  is  calculated  by  one 
of  two  formulas  given  in  Equations  11  or  12.  One  of  the  most  commonly 
used  formulations  is  by  Atkinson  and  Holliday  (1985), 

Um=  3.45(1010  -,pe)  (11) 

Knaff  and  Zehr  (2007)  evaluated  several  wind-pressure  relationships  and 
determined  that  Umax  calculated  using  Equation  11  had  a  negative  bias. 
Therefore,  the  functional  form  of  Koba  et  al.  (1990)  is  used  in  TWAVE  as 

AP  =  6.22  — 0.58(7 maJ.  -(£/„„/ 31.62)2  (12) 

The  user  can  either  use  Equations  n  or  12  by  selecting  the  toggle  centrpres 
in  the  TWAVE  Control  File  (Chapter  5). 

The  radius  of  maximum  wind  speed  Rm ax  is  used  to  estimate  the  size  of  the 
storm.  If  Rmax  is  not  available  (denoted  by  -99),  the  empirical  relation  of 
Willoughby  et  al.  (2006)  is  used  to  estimate  this  parameter  as 
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^max  =  46.4exp(-0.0155£/max+0.0169<p)  (13) 

where  tp  is  the  latitude  in  degrees.  Equation  13  is  valid  only  for  the 
Northern  Hemisphere.  Storm  track  positions,  wind  speeds,  central 
pressures,  and  radius  of  maximum  wind  speed  are  read  at  6-hr  intervals 
and  interpolated  hourly  estimates  are  calculated  by  a  cubic-spline  fit. 
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3  Astronomical  Tides 

The  astronomical  tides  affect  water  levels,  which  in  turn  influence  the 
storm  surges  and  waves  causing  flooding  of  island  shorelines.  In  TWAVE, 
tidal  water  elevations  are  computed  using  the  Oregon  State  Tidal 
Prediction  Software  OTPS  (http://www.coas.oregonstate.edu/research/ 
po/research/tide/otis.html),  which  operates  on  a  pre-calculated  tidal 
database.  Two  global  and  one  regional  (Hawaiian  Islands)  tidal  databases 
are  developed.  These  tidal  solutions  were  calculated  using  the  Oregon 
State  Inversion  Software  called  OTIS  (Egbert  and  Erofeeva  2002).  The  tide 
water  levels  and  velocities  are  calculated  with  a  least-squares  fit  of  along 
track  average  TOPEX/Poseidon  and  Jason  altimeter  data  to  the  Laplace 
Tidal  equations.  Eight  primary  (M2,  M4,  S2,  N2,  Ki,  Oi,  Pi,  and  Qi)  and 
two  minor  long -period  harmonic  constituents  (Mf  and  Mm)  are  provided 
as  complex  amplitudes  on  a  grid  with  V4-deg  resolution  for  the  global  tide 
database  (TPXO7.1)  and  1/30-deg  for  the  regional  Hawaiian  Islands  data¬ 
base  (Haw).  Additional  solutions  for  specific  regions  maybe  downloaded 
from  http://www.coas.oregonstate.edu/research/po/research/tide/ 

region. html,  and  placed  in  the  folder  called  OTIS/DATA/. 

The  OTPS  software  has  been  implemented  in  TWAVE  by  adding  programs 
that  prepare  OTPS  input  files  and  convert  output  files  named  Astronomi¬ 
cal  Tides  File  for  TWAVE  to  use.  These  ASCII  files  contain  hourly  tidal 
predictions  at  Nearshore  Stations.  Other  tidal  prediction  models  or  data¬ 
bases  may  be  incorporated  in  TWAVE  to  generate  the  Astronomical  Tide 
Files.  For  more  information  on  the  file  format,  how  to  specify  the  tidal 
solution  name,  and  other  input  parameters,  see  Chapter  5. 
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4  Waves 

Deepwater  wave  models 

Two  models  are  included  in  TWAVE  for  calculating  the  deepwater  waves. 
The  first  is  a  spectral  wave  model  and  the  second  is  a  parametric  model. 

Spectral  wave  model 

The  second-generation  spectral  wave  model  WISWAVE  (Wave  Informa¬ 
tion  Studies  WAVE  model)  is  implemented  in  TWAVE  to  simulate  deep¬ 
water  wave  generation,  growth,  propagation,  dissipation,  and  transforma¬ 
tion  (Hubertz  1992;  Resio  and  Perrie  1989).  WISWAVE  solves  the  discrete 
form  of  the  unsteady  wave  energy  balance  equation  in  two  dimensions  as 

f }  +  C,-™  =  tlS<  (14) 

where  N  is  the  wave  action  density  (energy  density/frequency),  t  is  time, 

Cg  is  the  group  velocity,  and  Si  are  the  source/sink  functions  for  wind 
input,  whitecapping  dissipation,  and  nonlinear  wave-wave  interactions. 
Tracy  and  Cialone  (2004)  showed  comparable  results  between  WISWAVE 
and  a  third-generation  wave  model  for  wave  hindcast  in  the  Gulf  of 
Mexico.  Wind  fields  are  input  at  l-hr  intervals.  The  model  output  at 
selected  grid  points  consists  of  wave  parameters  (wave  height,  peak  and 
mean  wave  period,  and  mean  wave  direction)  and  percent  total  energy  in 
each  frequency  band  of  the  wave  spectra.  The  grid  points  can  be  either 
entered  in  the  WISWAVE  Options  File  or  in  the  Offshore  Stations  File. 

The  Offshore  Stations  File  is  used  if  no  stations  are  specified  in  the 
WISWAVE  Options  File.  WISWAVE  is  setup  in  TWAVE  to  use  the  large 
PBL  model  grid. 

Parametric  wave  model 

The  WISWAVE  spectral  wave  model  is  the  most  computationally  intensive 
model  in  the  TWAVE  package.  A  simpler  and  faster  model  is  adequate  for 
planning  studies  and  for  simulation  of  hypothetical  storms.  Therefore,  a 
parametric  hurricane  wave  model  (Bretschneider  1972)  is  also  imple¬ 
mented  in  TWAVE.  This  model  assumes  a  constant  wave  height 
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distribution  based  on  a  slow  moving  hurricane  assumption.  The  deepwater 
maximum  wave  height  and  peak  period  are  calculated  every  hour  based  on 
the  following  empirical  relations  at  the  point  of  maximum  winds  as 


Hs  =5. 03  exp 


R  A P) 

max 

0.29a vf 

1  1  1 

4700  „ 

(15) 


8.6exp 


R  A P) 

max 

0.145am 

1  1  s 

9400  J 

(16) 


where  Ur  is  the  maximum  sustained  wind  speed  in  m/sec  at  10-m  eleva¬ 
tion,  Rmax  is  the  radius  of  maximum  wind  in  km,  A P  is  the  difference 
between  the  normal  atmospheric  pressure  and  the  central  pressure  of  the 
hurricane  specified  in  mm  of  mercury,  v /is  the  forward  velocity  of  the 
hurricane  in  m/sec  and  a  is  an  empirical  coefficient  assumed  to  be  l.o.  The 
spatial  distribution  of  wave  height  is  calculated  using  the  contour  values 
and  the  coordinates  with  respect  to  the  forward  direction  of  the  storm 
(Figure  2). 


Figure  2.  Normalized  deepwater  significant  wave  height  for  slow  moving  hurricane. 
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The  wave  period  is  estimated  as 

Ta  =12.1  jHJg  (17) 

The  wave  direction  is  perpendicular  to  the  radial  from  the  hurricane  center 
(Figure  2).  Although  the  parametric  wave  model  is  fast,  it  may  not  be 
appropriate  for  forecasting  or  hindcasting  storms  because  it  approximates 
many  physical  processes  governing  the  generation  and  growth  of  deep¬ 
water  waves.  It  does  not  account  for  wave  propagation  and  therefore 
underestimates  waves  that  are  more  than  approximately  3  Rmax  away  from 
the  storm  center.  It  is  recommended  that  this  model  only  be  used  for  sta¬ 
tions  near  the  storm  eye  and  for  hypothetical  events  for  which  the  uncer¬ 
tainty  of  input  parameters  is  greater  than  those  of  historical  events. 

Deep-  to  shallow-water  wave  transformation 

The  deepwater  wave  estimates  from  offshore  are  transformed  to  nearshore 
using  the  Wave  Information  Study  (WIS)  Phase-3  transformation  model 
WAVTRAN  (Jensen  1983;  Gravens  et  al.  1991).  The  wave  spectra  are  trans¬ 
formed  based  on  the  shoreline  orientation  with  respect  to  the  incident 
waves,  the  change  in  water  depth,  and  sheltering  effects.  WAVTRAN 
assumes  straight  and  parallel  bathymetric  contours  and  uses  an  input 
TMA  spectral  shape  (Bouws  et  al.  1985).  Sheltering  is  considered  by 
removing  the  wave  energy  from  sheltered  grid  points  (wave  directions  that 
are  sheltered  by  land).  The  directional  spreading  function  is  assumed  as  a 
cosine  to  the  4th  power  angular  distribution. 

Nearshore  wave  models 

The  term  nearshore  here  refers  to  the  coastal  region  from  intermediate 
(0.05  <  h/L  <  0.25)  to  shallow  water  ( h/L  <  0.05)  where  h  is  water  depth 
and  L  is  wavelength,  and  in  which  waves  shoal,  refract,  and  break.  The 
nearshore  water,  relative  to  the  mean  sea  level  (MSL),  includes  a  still  water 
level  (SWL)  component  (astronomical  and  barometric  tide  plus  wind- 
driven  surge)  and  wave  setup  component  (driven  by  momentum  lost  from 
the  waves  as  they  break.  Figure  3  shows  a  schematic  of  a  fringing  reef  and  a 
typical  wave  setup  distribution  across  the  reef.  Wave  runup  is  the  elevation 
of  highest  wave  excursion  relative  to  the  SWL  (thus,  it  includes  the  setup). 
Inundation  refers  to  the  horizontal  distance  landward  from  the  shoreline 
to  the  farthest  inland  reach  of  the  water  level  with  the  wave  runup. 
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Note  that  the  runup  and  inundation  are  defined  with  respect  to  the  SWL 
while  the  max  water  level  is  equal  to  the  runup  level  with  respect  to  the 
MSL. 

The  nearshore  wave  transformation  is  calculated  at  hourly  intervals  on 
predefined  transects  that  are  approximately  normal  to  shore.  There  are 
four  options  for  calculating  wave  heights  and  water  levels  over  the  reef.  All 
methods  use  the  same  incident  nearshore  wave  heights  and  water  levels, 
but  vary  in  their  complexity  and  the  amount  of  input  information 
required.  The  first  two  methods  are  used  when  the  bathymetric  profile  is 
not  available.  Method  l  calculates  the  wave  height  over  the  reef  by  assum¬ 
ing  a  stable  wave  height  and  calculates  reef  top  water  levels  using  simple 
empirical  formulas  for  wind  and  wave  setup.  The  second  method  calcu¬ 
lates  the  cross-shore  wave  height  variation  using  a  simple  energy  flux 
method,  but  also  uses  the  same  empirical  formulas  as  Method  l  to  calcu¬ 
late  reef  top  water  levels.  The  third  method  is  a  steady-state  model  which 
solves  the  lD  wave  energy  conservation  and  momentum  equations.  The 
last  method  is  a  time-dependant  nonlinear  wave  model  based  on  the 
Boussinesq  equations  (Demirbilek  et  al.  2008;  Demirbilek  and  Nwogu 
2007).  Methods  1-3  are  calculated  within  the  program  reefivave_TW.exe 
at  hourly  intervals  while  Method  4,  the  Boussinesq  wave  model,  is  run 
using  a  separate  batch  file  (run_boussid.bat).  These  methods  are  briefly 
discussed  next. 

Method  1:  Empirical  formulas 

Because  many  coral  reefs  are  characterized  by  wide  shallow  platforms, 
wave  heights  over  the  reef  are  usually  depth  limited.  Therefore,  the  wave 
height  estimate  will  also  be  fairly  accurate  as  long  as  the  mean  water  depth 
over  the  reef  top  can  be  estimated  reasonably  accurately.  The  setup  on  the 
reef  platform  is  separated  into  wave  and  wind  components 
Ar|  =  Aq  wave  +  Ar|  Wmd.  The  procedure  for  calculating  wave  setup  over  the 
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reef  is  as  follows.  First,  the  reef  platform  water  depth  is  adjusted  for  wave 
setup  using  an  empirical  equation  and  then  the  wind  setup  is  calculated 
using  semi-analytical  equation. 

Earlier  research  studies  (Seelig  1983;  Gourlay  1996a,  1996b)  have  shown 
that  the  reef  top  wave  setup  decreases  with  the  still  water  depth,  increases 
with  the  offshore  wave  power,  and  decreases  with  increasing  wave  steep¬ 
ness.  Using  these  findings,  Sanchez  et  al.  (2007)  developed  an  expression 
to  estimate  wave  setup  over  fringing  reef  profiles  similar  to  that  of  Seelig 
(1983)  as 


A'l„„  =  -0.18/!,.+0.481og10(//„2r)-5.53//o/r2  (18) 

where  hr  is  the  still-water  depth  on  the  reef  top.  For  fringing  reefs,  most  of 
the  wave  setup  occurs  in  a  narrow  region  near  the  reef  edge  where  the 
waves  break  and  lose  a  large  portion  of  their  energy  and  momentum.  The 
wind  setup  over  the  reef  platform  can  be  estimated  from  the  momentum 
balance  over  the  reef  as 


P  9dwr^r-  =  iwx  (19) 

ox 

where  dwr  =  hr  +  Ar\waVe  and  twx  is  the  wind  stress  in  the  x-direction 
defined  as  xwx  =  paCd\U\  U.  Here  pa  is  the  air  density,  U is  the  wind  velocity 
in  the  x-direction  (cross-shore),  and  Cd  is  the  wind  drag  coefficient.  There¬ 
fore,  by  assuming  that  Ar\waVe  remains  constant  over  the  reef  top  and 
following  Bretschneider  (1966),  Equation  19  can  be  integrated  over  the 
width  of  reef  W  to  give  an  expression  for  the  wind  setup  over  the  reef  plat¬ 
form  given  by 


ArUw  -  ~dwr  +  V asWU2  +  d2wr  (20) 

where  as  ~  2paCd/(gp).  Although  the  wind  setup  is  usually  small  compared 
to  the  wave  setup,  during  storms  and  for  wide  shallow  reefs,  the  wind 
setup  can  be  significant.  Figure  4  shows  a  comparison  of  measured  and 
calculated  setup  over  a  fringing  reef  profile  from  laboratory  experiments  of 
winds  and  waves  (Demirbilek  et  al.  2007).  In  this  reef  example,  the  wave 
setup  has  been  calculated  with  Equation  18  using  waves-only  experimental 
conditions.  The  wind  setup  from  Equation  20  is  added  to  calculate  the 
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total  setup  on  the  reef  top.  As  can  be  seen  from  Figure  4,  including  the 
wind  setup  significantly  improves  the  comparison  between  predicted  and 
measured  setup  on  the  reef.  Calculated  wave  setup  for  the  reef  profiles 
used  by  Seelig  (1983)  and  Gourlay  (1994)  are  also  included  in  Figure  4. 


Figure  4.  Comparison  of  various  empirical  formulas  for  setup  over  fringing  reef  to  laboratory 
measurements  of  Demirbilek  et  al.  (2007). 


The  stable  wave  height  over  the  reef  is  estimated  from  the  depth-limited 
y  =  Hrms/d  ~  0.4  (Dally  et  al.  1985).  This  estimate  may  not  be  valid  in 
many  applications  such  as  relatively  large  waves  over  narrow  reefs  or 
conversely  relatively  small  waves  over  wide  reefs.  The  width  of  the  reef 
platform  is  considered  in  Method  2.  Method  1  is  used  when  the  user 
specifies  a  water  depth  over  the  reef  top  but  does  not  provide  a  reef 
platform  width. 

Method  2:  Simplified  geometry  and  wave  model 

In  situations  where  the  reef  is  narrow  or  very  wide,  it  is  important  to 
consider  the  spatial  variation  of  wave  height  over  the  width  of  the  reef.  In 
this  approach,  the  wave  and  wind  setup  are  calculated  as  in  Method  1  to 
obtain  the  water  level  over  the  reef.  The  wave  height  is  then  propagated 
over  the  reef  using  a  simple  energy  flux  method  and  the  wave  breaking 
formulation  of  Dally  et  al.  (1985)  (see  description  of  Method  3  for  further 
details).  In  this  approach  the  water  depth  over  the  reef  is  not  coupled  to 
the  wave  energy  flux.  This  is  the  approach  used  in  Method  3.  The  wave 
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height  at  the  offshore  end  of  the  reef  top  can  be  either  set  to  a  breaking 
wave  height  using  the  empirical  equation  of  Komar  and  Gaughan  (1973) 
where  Hb  =  0.39 (H0/L0y°-2  (where  H0  is  deepwater  wave  height  and  L0  is 
deepwater  wave  length),  or  set  to  the  nearshore  wave  height.  This  option  is 
controlled  using  the  toggle  ibreak,  within  the  TWAVE  Control  File.  See 
Chapter  5  for  further  details  on  TWAVE  Control  File. 

Method  3:  Wave  energy  flux  model  (WAV1D) 

A  number  of  different  type  of  lD  models  based  on  the  conservation  of  wave 
energy  flux  and  momentum  have  been  used  extensively  and  successfully  in 
calculating  wave  transformation  and  water  levels  in  the  surf  zone 
(Gerritsen  1980;  Thornton  and  Guza  1983;  Dally  et  al.  1985;  Dally  1992; 
Larson  and  Kraus  1991;  Massel  and  Brinkman  2001;  Grasmeijer  and 
Ruessink  2003).  There  are  two  approaches  that  can  be  used  in  lD  energy 
flux  models.  Single  wave  models  (e.g.,  Dally  et  al.  1985)  simulate  the  wave 
transformation  of  a  representative  or  characteristic  wave  height,  wave 
period,  and  direction,  while  probabilistic  models  simulate  wave  transfor¬ 
mation  based  on  a  discrete  number  of  wave  classes  (e.g.,  Larson  and  Kraus 
1991).  Grasmeijer  and  Ruessink  (2003)  compared  single  wave  and  proba¬ 
bilistic  models  to  laboratory  and  field  measurements.  They  concluded  that 
the  single  wave  models  was  simpler  and  faster,  and  produced  similar  error 
statistics  to  the  probabilistic  models.  A  single  wave  approach  is  used  in 
TWAVE. 


Assuming  alongshore  uniformity  in  the  bathymetry,  the  time-averaged 
wave  energy  balance  equation  may  be  written  as 


~C^EC  cosG)  =  Dh  +  Df 
ox 


(21) 


where  E  is  the  wave  energy  density  per  unit  area,  Cg  is  the  wave  group 
velocity,  0  is  the  incident  wave  angle  relative  to  shore  normal,  and  Db  and 
D/are  the  wave  energy  dissipation  due  to  breaking  and  bottom  friction, 
respectively.  By  assuming  that  the  individual  wave  heights  follow  a 
Rayleigh  distribution  (Headquarters  (HQ)  USACE  2006)  given  by 


2  H 

p{H)  =  < exp 


H 


H.. 


(22) 
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where  Hrms  is  the  root-mean-squared  (rms)  wave  height,  and  the  wave 
energy  density  is  then  defined  as 


.  oo  . 

E  =  -pgf  H2p(H)dH  =-p  gH2ml 


(23) 


The  significant  wave  height  is  estimated  as  Hs  =  V2 Hrms .  The  wave  dissi¬ 
pation  due  to  bottom  friction  is  calculated  as 


fi fw 

2  nHrms 

1231 

Tsinh(kd) 

(24) 


where  fw  is  the  wave-related  friction  coefficient  given  by 


/u,  =  min  °-3,  exp  -6.0 +  5.2 (Ab/kw) 


-0.19 


(25) 


where  Ab  is  the  near-bed  wave  orbital  excursion  calculated  from  linear 
wave  theory  and  kw  is  a  hydraulic  roughness  length  assumed  to  be  equal  to 
the  physical  roughness  (Nielson  1992).  Lowe  et  al.  (2005)  estimated  the 
hydraulic  roughness  along  several  transects  of  a  coral  reef  and  found 
kw  =  0.16  ±  0.03.  Gerritsen  (1980)  suggested  values  between  0.025-0.125. 
Demirbilek  and  Nwogu  (2007)  considered  reef  roughness  in  terms  of  the 
Chezy  friction  factor,  and  their  equivalent  values  to  thefw  friction  coeffi¬ 
cient  ranged  from  0.01  to  0.15. 

Because  the  wave  breaking  term  in  Equation  21  is  the  dominant  term  in 
the  surf  zone,  the  estimated  wave  results  are  sensitive  to  the  wave  breaking 
formulation  being  used.  Two  types  of  formulations  can  be  used  to  estimate 
the  breaking  wave  dissipation.  The  first  formulation  adopted  for  wave 
breaking  dissipation  is  based  on  the  work  of  Battjes  and  Janssen  (1978). 
Recently,  Alsina  and  Baldock  (2007)  and  Janssen  and  Battjes  (2007)  inde¬ 
pendently  obtained  the  following  expression  for  wave  breaking  dissipation 
as 


D,  =  ^-pgB-1  H,ms 


16 


T  d 


(26) 
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where  B  is  the  wave  breaking  intensity  factor,  and  the  fraction  of  broken 
waves  Qb  is  estimated  as 

Qh  =  1  +  -4=( R3  +  -f^|exp(-i?2 )  -  erf  (R)  (27) 

3y/7l{  2  J 

In  Equation  27,  erf  is  the  error  function  and  R  =  Hb/Hrms,  and  the  breaking 
wave  height  is  obtained  from  Hb  =  0.88  tanh  (ykd/o.88)/k.  The  breaker 
index  y  is  calculated  using  the  expression  of  Battjes  and  Stive  (1985) 
y  =  Hb/d  =  0.5  +  0.4  tanh  (33  H0/L0 )  (Alsina  and  Baldock  2007),  where  d 
is  local  water  depth,  and  k  is  wave  number. 

Equation  26  is  obtained  by  following  a  similar  approach  as  in  Battjes  and 
Janssen  (1978)  without  assuming  a  truncated  Rayleigh  distribution  (e.g., 
using  a  complete  Rayleigh  distribution).  The  default  coefficients  in 
Equation  26  may  be  changed  in  the  TWAVE  Control  File.  Alsina  and 
Baldock  (2007)  compared  Equation  26  to  the  breaking  formulations  of 
Thornton  and  Guza  (1983),  Baldock  et  al.  (1998),  and  Ruessink  et  al. 
(2003)  and  found  that  Equation  26  produced  small  errors  for  steep  slopes. 
Sanchez  et  al.  (2007)  compared  several  breaking  formulations  for  coral 
reef  profiles  and  found  that  Equation  26  produced  the  best  results.  There¬ 
fore,  Equation  26  is  the  default  breaking  formulation  implemented  in 
TWAVE.  The  second  breaking  formulation  used  in  TWAVE  is  the  empiri¬ 
cal  model  of  Dally  et  al.  (1985)  expressed  as 

I\  =  |[£C5-min(£,Ei,)C!,]  (28) 

where  k  is  an  empirical  decay  coefficient  and  Es  is  the  energy  associated 
with  the  stable  wave  height  defined  as 

E,  =  ip9(rd)2  (29) 

o 

where  T  is  the  stable  wave  height  normalized  by  the  water  depth.  Follow¬ 
ing  earlier  works  (Dally  et  al.  1985;  Smith  1993;  Demirbilek  and  Panchang 
1998;  Zhao  et  al.  2001),  we  adapt  0.15  and  0.42  for  the  values  of  k  and  F. 
This  approach  has  been  used  successfully  in  modeling  wave  transforma¬ 
tion  over  irregular  bathymetry  including  reefs  (Gerritsen  1980;  Dally 
1992). 
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Equation  28  does  not  depend  on  the  wave  period  and  height  distributions. 
The  wave  height  distribution  is  used  only  to  obtain  a  characteristic  wave 
height  from  the  wave  energy.  Equation  28  works  best  for  monochromatic 
waves  because  it  was  developed  from  laboratory  studies  of  monochromatic 
waves. 

Equation  21  is  solved  simultaneously  with  a  time-averaged  and  depth- 
integrated  lD  momentum  equation  to  obtain  the  combined  wind  and  wave 
setup  over  the  reef-top  and  the  model  is  called  WAViD  (Demirbilek  et  al. 
2008).  Assuming  uniformity  in  the  alongshore  direction  (y-direction)  and 
ignoring  mixing  and  cross-shore  currents,  the  lD  cross-shore  momentum 
equation  (Dean  and  Dalrymple  1991)  may  be  written  as 

pg(h  +  vi)^-  +  ^^  =  Twx  (30) 

ox  ox 

where  r\  is  the  wave  setup  relative  to  the  SWL,  h  is  the  still-water  depth 
relative  to  a  specific  vertical  datum,  and  Sxx  is  the  wave  radiation  stress  in 
the  cross-shore  direction.  The  wave  radiation  stress  is  approximated  using 
linear  wave  theory  for  an  arbitrary  wave  angle  (Longuet-Higgins  and 
Stewart  1964)  as 


Sxx  =  E(n(  cos2  0  + 1)  -  0.5)  (31) 

where  0  is  the  incident  wave  angle  relative  to  shore  normal  and 
n  =  0.5  [1  +  2kd/sinh(2kd)].  Equations  21  and  30  are  solved  simultane¬ 
ously  from  deep  to  shallow  water. 

Demirbilek  et  al.  (2008)  provide  a  comprehensive  description  of  the 
WAViD  model,  with  four  example  applications  that  illustrate  the  model’s 
performance  for  wave  shoaling,  refraction,  breaking,  and  bottom  friction 
compared  to  data.  Sanchez  et  al.  (2007)  compared  the  energy  flux  model 
with  two  laboratory  experiments  of  wave  transformation  and  wave  setup 
over  fringing  reef-type  profiles.  They  reported  average  percent  errors 
(average  error  divided  by  offshore  wave  height)  of  less  than  25  percent  for 
wave  heights  and  setup  over  the  reef  platform  for  both  waves  cases  with 
strong  winds.  Although  WAViD  cannot  represent  nonlinear  transfer  of 
energy  in  the  wave  spectrum,  it  can  describe  the  average  wave  energy  dis¬ 
sipation  and  wave  setup  over  the  reefs.  WAViD  is  efficient,  robust,  and 
works  for  extreme  wave  conditions.  This  makes  WAViD  a  useful 
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engineering  tool  for  nearshore  wave  transformation  and  water  levels  along 
cross-shore  transects. 

Method  4:  Boussinesq  wave  model 

BOUSS-iD  is  a  time-dependent,  nonlinear,  phase-resolving  wave  model 
appropriate  for  estimating  nearshore  waves,  including  wave  transforma¬ 
tion  over  reefs,  wave  setup  and  wave  runup.  BOUSS-iD  is  a  lD  version  of 
the  Boussinesq  model  BOUSS-2D  (Nwogu  and  Demirbilek  2001;  Demir- 
bilek  et  al.  2005a  and  b;  Demirbilek  and  Nwogu  2007;  Demirbilek  et  al. 
2008).  The  model  solves  for  the  fully  nonlinear  Boussinesq  equations 
using  a  finite-difference  method  and  implements  a  one-equation  turbu¬ 
lence  closure  model  to  simulate  wave  breaking  in  the  surf  zone 
(Demirbilek  and  Nwogu  2007;  Demirbilek  et  al.  2008). 

BOUSS-iD  models  nearshore  wave  transformation  and  can  simulate  wave 
shoaling,  reflection,  bottom  friction,  nonlinear  wave-wave  and  wave- 
current  interactions,  wave  breaking,  wave  setup,  and  wave  runup  and 
overtopping  of  structures.  Within  TWAVE,  the  incident  waves  in 
BOUSS-iD  are  assumed  to  be  normal  to  shore  (wave  refraction  is  not 
considered).  The  wind  setup  is  included  by  adding  a  constant  water  level 
estimated  by  solving  Equation  21  across  the  reef  profile.  BOUSS-iD  is  the 
most  accurate  of  the  four  methods  to  calculate  nearshore  wave  transfor¬ 
mation  and  runup.  However,  it  is  also  the  most  computationally  intensive 
because  it  is  a  nonlinear  and  time-dependant  model.  BOUSS-iD  inputs  are 
a  bathymetry  file  (*.xy)  and  a  script  file  (*.gbat)  that  contains  the  model 
input  parameters,  including  wave  and  water  level,  and  the  post  processing 
commands. 

Because  users  of  TWAVE  will  most  likely  be  concerned  with  the  maximum 
coastal  inundation  and  not  its  time  history,  it  is  not  necessary  or  practical 
to  run  BOUSS-iD  at  hourly  intervals  for  each  transect  as  in  the  case  of  the 
previous  methods.  BOUSS-iD  is  run  only  for  the  time  periods  of  estimated 
maximum  coastal  inundation  (runup)  from  the  previous  methods.  See  the 
BOUSS-iD  Input  Files  section  in  Chapter  5  for  information  on  how  to 
prepare  BOUSS-iD  script  files.  BOUSS-iD  is  a  computationally  efficient 
model,  and  a  typical  run  time  for  an  8oo-m-long  transect  is  approximately 
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Calculation  of  wave  runup  statistics 

Two  methods  are  used  in  the  estimation  of  wave  runup  and  inundation 
statistics  in  TWAVE.  The  first  method  uses  empirical  formulas  embedded 
within  the  WAViD  model.  The  second  method  applies  a  zero-crossing 
analysis  to  BOUSS-iD  runup  time  series. 


Empirical  formulas 

For  design  and  planning  purposes,  it  is  useful  to  describe  the  runup  in 
terms  of  the  maximum  Rmax  and  2-percent  exceedance  level  R2%  (runup 
exceeded  by  2  percent  of  waves).  R2%  is  estimated  in  three  different  ways 
The  first  uses  an  empirical  formula  from  HQUSACE  (2006), 

Equation  VI-5-7,  applicable  for  irregular  waves  on  a  composite  slope  as 


R2%  _  [1-5^YrYbYftYp  for  0.5  <2 

Hs  [3.0YrYhYp  for  ^  >  2 


(32) 


where: 

y r  =  surface  roughness  factor 

Y b  =  berm  influence  factor  (=  1) 

jh  =  wave  height  distribution  factor  (=1) 

Yp  =  incident  wave  angle  factor 

c,  =  surf  similarity  parameter  given  by  £  =  tana /  ^HJ  L0  where 
Hr  is  the  wave  height  over  the  reef  top. 

The  surface  roughness  factor  is  approximately  equal  to  1.0  for  sandy 
beaches.  For  more  information  on  the  roughness  factor  refer  to  the 
Coastal  Engineering  Manual  Table  VI-5-3  (HQUSACE  2006).  On  fringing 
reefs,  waves  usually  break  near  or  on  the  reef  rim,  and  waves  may  reform 
again  over  the  reef  top.  This  suggests  that  wave  height  over  the  reef  top 
should  be  used  in  Equation  32  instead  of  the  deepwater  wave  height.  Wave 
transformation  over  reefs  is  characterized  by  the  presence  of  infragravity 
waves  generated  by  non-linear  wave  energy  transfer.  For  this  reason,  the 
deepwater  wavelength  is  used  instead  of  the  wavelength  over  the  reef.  The 
incident  wave  angle  factor  Yp  is  calculated  as 
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Sanchez  et  al.  (2007)  found  that  R2%  and  Rmax,  calculated  with  the  empiri¬ 
cal  runup  equations  by  Mase  (1989),  were  within  a  50  percent  error  of 
laboratory  measurements  for  fringing  profile.  Although  these  empirical 
methods  are  approximate,  they  provide  a  qualitative  and  robust  means  of 
estimating  the  largest  coastal  inundation.  The  time  of  maximum  coastal 
inundation  or  runup  may  not  necessarily  occur  during  the  period  of  largest 
coastal  waves,  due  to  the  combined  influence  of  storm  surge,  tides,  wave 
period  and  direction,  and  therefore  may  vary  for  different  locations  around 
an  island.  Once  the  critical  time  period  is  identified  for  each  transect,  the 
user  may  choose  to  run  BOUSS-iD  for  these  times.  The  TWAVE  variable 
ir2per  defines  which  empirical  equations  for  R2%  shall  be  used.  Options 
are:  o  for  Mase  (1989);  1  for  HQUSACE  (2006);  2  for  Hedges  and  Mase 
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(2004).  Similarly,  irmax  defines  the  equations  for  Rnmx  with  options:  o  for 
Mase  (1989),  1  for  Seelig  (1983).  See  Chapter  5  for  more  details  on  how  to 
use  input  variables. 

Analysis  of  wave  runup  time  series 

In  the  second  method,  the  2-percent  exceedance  level  runup  R2%  and 
maximum  runup  Rmax  for  each  transect  are  calculated  from  BOUSS-iD 
runup  time  series.  This  analysis  consists  of  an  upward  zero-crossing 
analysis  to  identify  individual  runup  events.  Noise  near  the  SWL  is  elimi¬ 
nated  by  setting  a  minimum  amplitude  and  period  for  runup  events  equal 
to  5  and  25  percent  of  the  nearshore  wave  height  and  period,  respectively. 
Runup  values  are  sorted  in  ascending  order  and  percentiles  are  estimated 
as  Pi  =  100  (i  +  0.5 )/n,  where  Pi  is  the  i-th  percentile  and  n  is  the  number 
of  runup  events.  The  98th  percentile  is  linearly  interpolated  to  obtain  the 
2-percent  exceedance  level  runup.  The  maximum  runup  is  calculated 
directly  from  the  time  series.  Inundation  statistics  are  calculated  by  deter¬ 
mining  the  location  of  maximum  water  level  and  runup  along  the  transect 
elevation  profile. 
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5  TWAVE  Modeling  System 

This  chapter  explains  how  to  use  the  TWAVE  package.  The  first  section 
describes  the  functionality  of  the  Excel  interface  developed  for  TWAVE. 
The  subsequent  sections  provide  a  detailed  description  of  the  input  and 
output  files. 

Excel  interface 

TWAVE  modeling  system  is  implemented  in  Microsoft  Excel®,  which 
serves  as  the  graphical  user  interface  (GUI)  for  the  system’s  visualization 
of  input  and  output  files.  Excel  was  chosen  as  an  interface  because  it  is 
available  on  most  personal  computers  (PCs)  and  avoids  using  other  pro¬ 
prietary  software.  The  users  of  TWAVE  do  not  need  any  commercially 
available  compilers  or  plotting  programs,  and  once  the  package  is 
unzipped,  it  is  ready  to  be  used.  Separate  Excel  files  are  used  for  individual 
island  applications  for  easier  data  management  and  organization.  Each 
Excel  files  contains  a  number  of  worksheets  for  different  functions.  The 
worksheets  are  color-coded  for  easy  navigation.  Figure  5  shows  the  first 
worksheet  of  the  TWAVE  package  which  contains  the  names  and 
descriptions  of  worksheets. 

Users  of  TWAVE  should  familiarize  themselves  with  the  Table  of  Contents 
and  each  of  the  worksheets  before  attempting  to  use  it  in  practical  applica¬ 
tions.  The  worksheets  are  arranged  sequentially  in  the  order  they  are 
utilized  and  the  tab  colors  are  coded  based  on  the  worksheet  functions. 
Yellow  tabs  correspond  to  worksheets  with  reference  information  such  as 
the  table  of  contents  and  maps.  Green  tabs  represent  worksheets  used  to 
visualize  input.  Red  tabs  are  used  for  worksheets  that  contain  output 
results  in  tables  and  plots.  The  black  tab  named  STR  is  the  steering  work¬ 
sheet  that  contains  the  hyperlinks  to  run  specific  TWAVE  applications 
(Figure  6).  The  applications  include  pre-  and  post-processing  routines, 
numerical  models,  and  batch  files  (*.bat).  Batch  files  are  scripts  that  run 
several  commands  and/or  executables. 
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Figure  5.  Screen  snap  shot  of  first  worksheet  in  TWAVE. 


Before  executing  any  models  available  in  the  TWAVE,  the  current  path  in 
Excel  must  be  set  to  the  TWAVE  work  directory.  This  is  done  by  using  the 
Open  dialog  in  Excel  which  can  be  accessed  by  clicking  on  File  \  Open  or 
the  using  the  Open  button  The  current  path  is  set  by  browsing  to  the 
directory  /TWAVE/ work  and  then  clicking  on  Cancel  once  the  correct  path 
is  in  the  dialog  box.  If  this  procedure  is  repeated,  the  default  path  on  the 
Open  dialog  box  would  be  set  to  /TWAVE/work. 

To  execute  an  individual  program  in  the  TWAVE,  click  on  the  respective 
box  in  the  Steering  worksheet  STR  (Figure  6).  This  should  pop  up  a  DOS 
window  showing  the  screen  output  for  that  program  and  indicating 
whether  or  not  the  program  ran  successfully.  Once  the  program  is 
finished,  the  user  should  review  the  screen  information,  check  for  any 
error  messages,  and  press  any  key  to  close  the  window. 
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Figure  6.  Screen  snap  shot  of  steering  worksheet  (each  white  box  represents  a  hyperlink). 

After  running  a  specific  model,  the  results  can  be  viewed  in  Excel  by 
updating  the  external  data  to  be  read  into  specific  worksheets  from  output 
files.  This  is  done  by  right-clicking  on  a  cell  with  external  data  and  clicking 
on  the  menu  Data  \  Refresh  Data  or  clicking  on  the  Refresh  button  I  . 
Refreshing  data  will  also  automatically  update  the  associated  plots. 
Worksheets  used  to  view  output  files  can  be  used  as  templates  for  viewing 
several  files  of  the  same  type  by  copying  a  worksheet  and  then  replacing 
the  external  data  file  names.  A  description  of  the  TWAVE  programs  and 
their  input  and  output  files  is  provided  in  Table  l.  Some  programs  such  as 
the  wiswave  TW.exe  are  run  with  the  batch  files. 


Table  1.  TWAVE  programs  and  associated  input  and  output  files. 


Program 

Description 

Input 

Output 

phm_TW.exe 

Parametric  Hurricane  Model 

TWAVE  Control  File 

Storm  Track  File 

Nearshore  Wind  and  Storm  Surge  Files 

pwm_TW.exe 

Parametric  Wave  Model 

TWAVE  Control  File 

Storm  Track  File 

Offshore  Waves  Files 

pblinput_TW.exe 

Prepares  the  input  files  for  the  PBL  model 

Run  within  the  batch  file  run_PBL.bat 

TWAVE  Control  File 

Storm  Track  File 

6-hr  Storm  Parameters  File 

1-hr  Storm  Parameters  File 

pbl_TW.exe 

Planetary  Boundary  Layer  (PBL)  model 
as  described  in  Thompson  and  Cardone  (1996) 
configured  on  a  WISWAVE  grid 

TWAVE  Control  File 

6-hr  Storm  Parameters  File 

1-hr  Storm  Parameters  File 

1-hr  Wind  and  Pressure  Fields 

wiswave_TW.exe 

Spectral  wave  growth  and  propagation  model  as 
described  in  Hubertz  (1992).  The  model  is  run 
within  the  batch  file  run_WISWAVE.bat 

TWAVE  Control  File 

1-hr  Wind  and  Pressure  Fields 
WISWAVE  Options  File 

Offshore  Station  Spectra  File 

Last  Significant  Wave  Height  File 

Maximum  Significant  Wave  Height  File 

spec_TW.exe 

Extracts  time  series  of  wind  and  wave 
parameters  for  each  WISWAVE  observation 
station  (offshore  station).  Run  within  the  batch 
file  run_WISWAVE.bat 

TWAVE  Control  File 

Offshore  Station  Spectra  File 
WISWAVE  Options  File 

Offshore  Station  Time  Series  Files 

wavetran_TW.exe 

Transforms  hourly  deepwater  wave 
parameters  to  shallow  water  to  represent 
near-breaking  waves  approaching  the  coast 
(Jensen  1983;  Gravens  et  al.  1991) 

TWAVE  Control  File 

Storm  Track  File 

Nearshore  Station  File 

Offshore  Station  File 

Offshore  Station  Time  Series  Files 

Nearshore  Wave  Files 

wi  nd_i  n  v_ba  r_TW.  exe 

Extracts  the  wind  and  atmospheric  pressure  from 
the  PBL  model  output  data  and  computes  the 
storm  surge  time  series  at  nearshore  stations 

TWAVE  Control  File 

Storm  Track  File 

Nearshore  Station  Information  File 

1-hr  Wind  and  Pressure  Fields 

Nearshore  Wind  and  Storm  Surge  Files 

(Continued) 

Key  to  filename  components:  ww  =  wis  station  #,  ssssss  =  storm  #,  xx  =  nearshore  station,  ppp  =  profile  #,  yyyy  =  year 

GO 

O 
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Table  1.  Concluded. 


Program 

Description 

Input 

Output 

prep_tide.exe 

Prepares  the  input  files  for  the  TPXO  tide  model 

TWAVE  Control  File 

Storm  Track  File 

Station  Information  File 

Tide  Model  Setup  File 

Tide  Model  Lat/Lon/Time  File 

OTPS.exe 

Computes  time  series  of  astronomical  tides  at 
nearshore  stations  using  the  Oregon  State  Tidal 

Prediction  Software  (Egbert  and  Erofeeva  2002) 

Tide  Model  Setup  File 

Tide  Model  Lat/Lon/Time  File 

OTPS  Tides  File 

proc_tide.exe 

Converts  the  OTPS  output  file  to  TWAVE  format 

TWAVE  Control  File 

OTPS  Tides  File 

Astronomical  Tides  File 

combine_TW.exe 

Prepares  the  input  file  for  setup_TW.exe  by 
combining  the  time  series  of  tides,  storm  surge 
and  wave  conditions  at  nearshore  stations 

TWAVE  Control  File 

Nearshore  Station  Information  File 

Tide  Model  Setup  File 

Nearshore  Wave  Files 

Nearshore  Wind  and  Storm  Surge  Files 
Nearshore  Tide  File 

Nearshore  Time-series  Files 

WAVlD.exe 

Computes  time  series  of  wave  levels,  wave 
ponding,  wave  setup,  runup  and  inundation 
at  transects  along  the  coast 

TWAVE  Control  File 

Nearshore  Station  Information  File 

Profile  Information  File 

Nearshore  Time-series  File 

Transect  Bathymetry  Files 

Transect  Time-series  Files 

Transect  Profile  Files 

WAV1D  Runup  Statistics  File 
run_boussld.bat 

boussld_v3plexe 

ID  Boussinesq  wave  model 

Run  within  the  batch  file  run_boussld.bat 

B0USS-1D  script  files 

Transect  Bathymetry  Files 

B0USS-1D  Transect  Files 
B0USS-1D  Runup  time  series 

proc_ts.exe 

Calculates  the  runup  statistics  from  B0USS-1D  results 

B0USS-1D  Runup  time  series 

B0USS-1D  Runup  Statistics 

File 

Key  to  filename  components:  ww  =  WIS  station  #,  ssssss  =  storm  #,  xx  =  nearshore  station,  ppp  =  profile  #,  yyyy  =  year. 

co 
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Input  files 

Input  parameters  for  different  models  used  inside  TWAVE  must  be  speci¬ 
fied  in  separate  ASCII  files.  Input  files  include  the  storm  track,  maximum 
sustained  wind  speed,  location  of  nearshore  stations,  and  cross-shore 
bathymetry  profiles.  Each  input  file  name  is  associated  with  a  variable. 
Table  2  contains  the  TWAVE  input  files  names,  file  description,  file  loca¬ 
tion,  and  corresponding  variable  name.  Input  file  names  must  be  specified 
within  the  Control  File  next  to  the  corresponding  variable  names.  An 
example  of  the  Control  File  is  provided  in  the  next  section.  All  input  data 
to  TWAVE  are  in  SI  units,  except  the  wind  speed  must  be  specified  in 
knots. 


Table  2.  TWAVE  default  input  files. 


File(s) 

Description 

Name 

Folder 

TWAVE 

Variable 

Control  File 

Files  used  to  input  file  names 
and  parameters  in  TWAVE 

twave.ctl 

input 

NA 

Best  Track  File 

Best  Track  File  issued  by  the 
Joint  Typhoon  Warning  Center 

NA 

best_track 

fpred 

Storm  Track  File 

Input  file  for  the  PBL  model 
which  contains  the  storm  track 
and  maximum  wind  speed 

track_ssssss_N  AM  E.txt 

work 

ftrack 

Offshore  Stations 
File 

Contains  the  location  and 
depth  of  the  offshore  stations 

offstajnfo.inp 

input 

foffsta 

Nearshore  Stations 
File 

Contains  the  coordinates  of  the 
nearshore  stations  and  sets 
various  parameters  related  to 
them 

nrstajnfo.inp 

input 

fnrsta 

Transects  File 

Contains  the  location  and 
orientation  of  nearshore 
bathymetry  profiles  and  sets 
various  parameters  related  to 
them 

transects.inp 

input 

ftrans 

WISWAVE  Options 
File 

Contains  the  settings  for  the 
WISWAVE  model 

options_wiswave.dat 

input 

fwisopt 

B0USS-1D  Script 
Files 

Contain  B0USS-1D  setup 
information  and  post¬ 
processing  commands 

Pppp_ssssss.gbat 

BousslD 

NA 

ssssss=  storm  ID,  NAME  =  storm  name  (optional),  ppp  =  transect  ID  number,  NA  =  not  applicable. 
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Control  File 

The  Control  File  defines  all  input  variables  for  TWAVE.  It  includes  file 
names  and  empirical  constants  for  various  numerical  models  imbedded  in 
TWAVE.  The  name  of  the  Control  File  is  twave.ctl  and  this  name  cannot  be 
changed.  An  example  twave.ctl  file  is  shown  below. 


TWAVE  Control  File 
Version  1 . 0 


START_CTL 


!For  Kauai 

ftrack  track_199210.txt 
foffsta  of f sta_kauai . inp 
fnrsta  nrsta_kauai . inp 
fprof  prof iles_kauai . inp 

fwisopt  options_kauai.dat 
region  hawaii 
fwind  wind.dat 


! Storm  Track  File 
! Offshore  Station  Information 
INearshore  Station  Information 
IProfile  Information  File 
1WISWAVE  options  file 
!  specifies  which  region  to  run 
!Wind  File 


windmodel  0  !0-Holland,  3-Rankine 

tidemodel  haw  !OTIS  tidal  solution 


File 

File 


END_CTL 


The  list  of  the  possible  variables  that  may  be  used  in  the  Control  File  is 
provided  in  the  Excel  worksheet  labeled  as  VAR.  Each  variable  must  be 
entered  in  a  separate  line  within  the  section  between  the  START_CTL  and 
END_CTL  by  specifying  the  variable  name  followed  by  its  value  or  file 
name.  Metadata  information  (e.g.,  version,  project  name,  etc.)  may  be 
specified  anywhere  by  placing  an  exclamation  sign  in  front  of  the  text. 
Table  3  shows  the  complete  list  of  variables  that  may  be  specified  in  the 
Control  File.  Table  4  provides  different  run  options  available  in  TWAVE. 

The  Control  File  may  be  directly  edited  with  Excel  by  clicking  on  its 
hyperlink  in  the  steering  worksheet  STR.  Use  a  default  program  to  open 
files  with  the  same  extension  for  the  hyperlink  to  work  properly.  If  the  file 
does  not  open  when  clicking  on  the  hyperlink,  go  to  Windows  Explorer, 
right-click  on  the  file,  click  on  properties  and  then  select  a  text  editor  as 
the  default  program  (e.g.,  Wordpad). 
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Table  3.  List  of  variables  for  TWAVE  Control  File. 


Name 

Description 

Default  File  Name  or 
Value 

ftrack 

Storm  Track  File  name 

NA 

fpred 

Name  of  Best  Track  File 

NA 

asave 

Date  and  time  at  which  a  warm  start  file  is 
saved  as  yyyymmddhh 

none 

aread 

Date  and  time  at  which  a  warm  start  file  is  read 
as  yyyymmddhh 

none 

fsnaps 

Name  of  intermediate  storm  data  file  at  6-hr 
intervals 

snaps.data 

fhisto 

Name  of  intermediate  storm  data  file  at  1-hr 
intervals 

histo.data 

fwind 

Name  of  storm  wind  and  pressure  file 

wind.dat 

fwmax 

Maximum  wind  field 

Wmax.dat 

fwisopt 

Name  of  WISWAVE  options  file 

options_wiswave.dat 

foffsta 

Name  of  Offshore  Station  Information  File 

offstajnfo.inp 

fnrsta 

Name  of  Nearshore  Station  Information  File 

nrstajnfo.inp 

fprof 

Name  of  Profile  Information  File 

profiles.inp 

fhsmax 

Maximum  significant  wave  height  at  every  grid 
point 

Hs_max.dat 

fhslast 

Significant  wave  height  at  every  grid  point  at  last 
computation  time  step 

Hs_last.dat 

ftidesetup 

Tide  Model  Setup  File 

tides_setup.inp 

ftide 

Name  of  file  containing  a  time  series  of  the  tide 

tide_otis.dat 

flatlontime 

Tide  Model  Lat/Lon/Time  File 

latlontime.txt 

pa 

Atmospheric  pressure  in  mbar 

1010.0 

Ka 

Conversion  factor  between  1-min  to  10-min 
winds 

0.87 

Km 

Conversion  factor  between  PHM  winds  and 

10-m  reference  height  winds 

0.8 

reefbeg 

Threshold  depth  above  which  a  reef  can  appear 
(feet,  MSL) 

-2 

gammab 

Ratio  of  Hs  to  water  depth  over  the  reef 

0.6 

gamma 

Non-dimensional  stable  RMS  wave  height 

0.4 

dx 

Spatial  resolution  of  WAV1D 

0.5 

fric 

Runup  reduction  factor  due  to  non-smooth  slope 

0.6 

B 

Wave  breaking  intensity  factor  for  wave  breaking 

1 

xRankine 

Shape  factor  for  modified  Rankine  vortex  wind 
model 

0.5 

ksr 

Nikuradse  roughness  factor  for  reef 

0.16 

kappa 

Empirical  wave  decay  coefficient  for  wave 
breaking 

0.15 
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Table  4.  TWAVE  run  options. 


Name 

Description 

Possible  Values* 

Default 

region 

Specifies  which  region  to  run 

Hawaii  or  Guam 

Hawaii 

windmodel 

Toggle  for  selecting  a  parametric 
wind  model 

0  -  Holland  model 

3  -  Modified  Rankine  vortex 

0 

centrpres 

Specifies  formulation  for  central 
pressure  when  not  given 

0  -  Koba  et  al.  1990 

1  -  Atkinson  and  Holliday 

1985 

0 

tidemodel 

Toggle  for  selecting  OTPS  tidal 
solution 

haw  -  Hawaiian  Islands 

tpxo6.2  -  Global  solution  6.2 
tpxo7.0  -  Global  solution  7.0 

tpxo6.2 

Bconst 

Toggle  for  specifying  a  constant 

B  factor  for  Holland  Model 

<0  -  Calculated 

>0  -  Constant  value 

-99 

Rmaxconst 

Toggle  for  specifying  a  constant 

Rmax 

<0  -  Calculated  (if  Rmax< 0) 

>0  -  Constant  value 

-99 

phm4wis 

Toggle  for  parametric  wind 
models  to  output  WISWAVE  wind 
field* 

1  -  Yes 

0-  No 

1 

phm4adcirc 

Toggle  for  parametric  wind 
models  to  output  an  ADCIRC 
compatible  wind  field 

1- Yes 

0-  No 

1 

ibreak 

Selects  the  method  for 
estimating  the  breaking  wave 
height 

0  -  Komar  and  Gaughan  1973 

1  -  Breaking  wave  height  equal 
nearshore  wave  height 

0 

iwave_diss 

Selects  the  method  to  calculate 
wave  dissipation  due  to  breaking 

0  -  Dally  et  al.  1985 

1  -  Alsina  and  Baldock  2007 

0 

iwave_fric 

Specifies  whether  to  use  wave 
bottom  friction 

0-No,  1-Yes 

1 

ir2per 

Selects  the  method  for 
calculating  the  2%  exceedance 
runup 

0  -  Mase  1989 

1-CERC 

2  -  Hedges  and  Mase  2004 

0 

irmax 

Selects  the  method  for 
calculating  the  maximum  runup 

0  -  Mase  1989 

1  -  Seelig  1983 

0 

*PBL  model  always  outputs  WISWAVE  wind  field. 

Storm  Track  File 

The  Storm  Ti'ack  File  in  TWAVE  provides  storm  track  positions,  wind 
speeds,  central  pressure,  and  radius  of  maximum  wind  speed  at  6-hr 
intervals.  The  file  name  is  “track_ssssss_ATAME.txt”,  where  ssssss  is  a 
storm  identifier  number  and  NAME  refers  to  the  name  of  the  tropical 
storm  or  hurricane.  The  name  of  the  storm  is  optional  but  the  six-digit 
storm  ID  is  required.  The  ID  convention  used  in  TWAVE  is  the  year 
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followed  by  the  storm  event  number.  This  allows  users  to  sort  storms  in  a 
chronological  order  in  folders  for  easier  data  management.  The  Storm 
Track  File  must  be  located  in  the  work  subdirectory.  The  name  of  the 
Storm  Track  File  is  specified  in  the  Control  File  using  the  variable  name 
ftrack.  An  example  of  the  format  of  a  Storm  Track  File  is  provided  below 
for  Hurricane  Iniki  (track_199210_INIKI.txt): 


yyyymmddhh  LatN  LonE 
1992090918 
1992091000 
1992091006 
1992091012 
1992091018 
1992091100 
1992091106 
1992091112 
1992091118 
1992091200 
1992091206 
1992091212 


W_knt  Pres_mb  Rmax_km 
14.116 
14 . 328 
14 . 532 
14 . 685 
14 . 918 
15 .217 
15.658 
16.368 
17 .078 
18.285 
19.765 
21.273 


13. 

.  4 

205. 

.  7 

80 

984 

13. 

.  8 

204  . 

.  5 

85 

980 

14  . 

.  3 

203. 

.  1 

90 

960 

14  . 

.  7 

202  . 

.2 

100 

960 

15. 

.2 

201 . 

.  4 

100 

951 

15. 

.  9 

200  . 

,  7 

110 

948 

16. 

.  8 

200  . 

.2 

115 

947 

18  . 

.2 

200  . 

.  0 

120 

939 

19. 

.  5 

200  . 

.  0 

125 

938 

21 . 

.  9 

200  . 

.29 

115 

945 

23. 

.  7 

200  . 

.  6 

100 

959 

25. 

.  7 

201 . 

.  0 

80 

980 

The  Storm  Track  File  must  have  one  header  line  followed  by  lines  consist¬ 
ing  of  six  columns  of  data  separated  by  blank  spaces.  The  header  line  has 
the  names  of  input  variables,  with  the  first  column  containing  the  date  and 
time  (GMT)  specified  as  yyyymmddhh,  where  yyyy  is  the  year,  mm  is  the 
month,  dd  is  the  day,  and  hh  is  the  hour.  The  latitude  and  longitude  values 
are  specified  in  the  second  and  third  columns  in  degrees  N  (north)  and  E 
(east),  respectively.  The  forth  to  sixth  columns  are  the  maximum  sustained 
wind  speed  in  knots,  central  pressure  in  mbar,  and  radius  of  maximum 
wind  speed  in  km,  respectively.  If  the  storm’s  central  pressure  or  the 
radius  of  maximum  wind  speed  is  not  known,  a  negative  value  should  be 
placed  in  the  respective  columns,  and  TWAVE  will  calculate  estimated 
values  of  these  parameters.  The  Storm  Track  File  may  be  viewed  in  Excel 
using  the  worksheet  TRK. 

Storm  databases 

Two  storm  databases  are  used  in  TWAVE.  The  first  storm  database  is  the 
NOAA  tropical  storm  and  hurricane  database  HURD  AT  (Landsea  2004) 
available  at  (http://www.aoml.noaa.gov/hrd/hurdat/).  This  storm  data¬ 
base  covers  the  Atlantic  Ocean,  Gulf  of  Mexico,  and  Caribbean  Sea  for  the 
time  period  from  June  1851  to  December  2007.  The  second  database  from 
the  Joint  Typhoon  Warning  Center  (JTWC),  Naval  Pacific  Meteorology 
and  Oceanography  Center,  is  the  western  North  Pacific  Ocean  database 
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fhttps: //metocph.nmci.naw.mil/jtwc/best  tracks/).  The  HURDAT  data¬ 
base  is  provided  in  TWAVE  format  as  Storm  Track  Files  located  in  the 
TWAVE  subdirectory  data\HURDAT.  The  JTWC  database  is  provided  in 
its  original  format  as  best  track  files  located  in  the  TWAVE  subdirectory 
data\JTWC.  The  JTWC  best  track  files  can  be  converted  to  TWAVE  Storm 
Track  Files  by  using  the  prep_besttrack_TW.exe  utility  available  in  the 
steering  worksheet.  The  Best  Track  files  should  be  placed  in  the  subfolder 
named  best_track,  and  their  names  specified  in  the  Control  File  using  the 
variable  nam  efpred. 

WISWAVE  Options  File 

This  file  contains  the  input  settings  for  the  WISWAVE  model.  The  file 
includes  WISWAVE  setup  information,  including  frequency  bins,  number 
of  directional  bins,  and  information  about  grid  cells  that  are  land  or  water. 
The  file  should  be  located  in  the  input  folder.  The  default  file  name  is 
options_wiswave.dat  and  may  be  changed  using  the  TWAVE  variable 
fwisopt.  For  additional  information  about  this  file,  the  user  is  referred  to 
the  WISWAVE  manual  (Hubertz  1992). 

Offshore  Stations  File 

The  wave  estimates  from  the  deepwater  wave  models  are  saved  at  the  user- 
specified  Offshore  Stations  File.  These  wave  conditions  can  be  trans¬ 
formed  to  shallow  water  for  nearshore  wave  estimates.  For  additional 
information,  see  the  section  Deep-  to  Shallow- Water  Wave  Transforma¬ 
tion.  The  locations  and  depths  of  Offshore  Stations  are  specified  in  the 
Offshore  Stations  File.  This  file  must  be  located  in  the  input  folder  and  the 
name  of  the  file  is  specified  in  the  Control  File  using  the  variable/ojf/sta. 

An  example  of  the  default  Offshore  Stations  File  for  the  island  of  Kauai,  HI 
( offsta_kauai.inp ),  is  shown  below: 


Of f Sta 

G01 

G02 

G03 

GO  4 

G05 

GO  6 

G07 

G08 

GO  9 

G10 


LatOf f 
22 . 083 
22 . 000 

21 . 917 

21 . 917 

21 . 917 

21 . 833 

21 . 833 

21 . 833 

21 . 833 

21 . 833 


LonOf f 

200 .167 

200 .167 

200 .167 
200.250 
200 . 333 
200 . 333 
200 . 417 
200 . 500 
200.583 
200 . 667 


DepOf f 
-999 
-999 
-999 
-999 
-999 
-999 
-999 
-999 
-999 
-999 
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The  wave  estimates  at  user-specified  Offshore  Stations  File  are  inter¬ 
polated  from  the  WISWAVE  grid  solution  file  using  the  nearest  neighbor 
interpolation  method.  The  offshore  depths  should  be  indicated  as  positive 
depths  or  as  -999  for  default  deepwater  stations.  The  Offshore  Stations 
File  may  be  viewed  in  Excel  using  the  worksheet  named  OFFSTA. 

Nearshore  Stations  File 

Hourly  estimates  of  storm  surge,  tides,  wind,  wave  heights,  wave  periods, 
and  wave  directions  (wave  conditions  optional)  are  output  at  user- 
specified  nearshore  stations.  These  stations  can  also  be  used  to  output 
water  levels  at  buoy  locations  located  in  shallow  water  or  stations  within 
harbors  where  water  elevations  are  desired  but  wave  heights  are  not.  In 
the  Nearshore  Stations  File,  the  user  may  specify  an  offshore  station 
Offsta  or  N/A  for  “Not  Available”  indicating  nearshore  wave  transforma¬ 
tion  is  not  necessary.  In  the  first  case,  the  waves  will  be  transformed  from 
the  specified  offshore  station  to  shallow  water  using  WAVETRAN  (Jensen 
1983;  Gravens  et  al.  1991). 

If  no  offshore  stations  are  specified,  wave  heights  are  set  to  zero  and  the 
Nearshore  Stations  File  will  output  the  water  level  only.  The  N/A  should 
be  used  as  a  placeholder  for  no  offshore  stations.  The  file  should  be  saved 
in  the  input  folder.  For  visualization  purposes,  the  information  in  this  file 
may  be  copied  to  the  Excel  worksheet  labeled  as  NRSTA.  The  default  name 
of  the  input  Nearshore  Stations  File  ( nrsta__info.inp )  can  be  changed  in 
the  TWAVE  Control  File  using  the  variable  mane  fnrsta.  An  example  of 
this  file  is  shown  below  for  the  island  of  Kauai: 


NrSta 

LatN 

LonE 

DepNr 

OffSta 

Opt 

AngBath 

KSH 

KSH1 

KSH2 

130 

21 .8990 

200 . 4058 

3 

NA 

1 

5 

0 

0 

0 

103 

21 . 8855 

200 . 4396 

32 

G07 

1 

20 

0 

0 

0 

128 

21 . 8777 

200 . 4594 

35 

G08 

1 

20 

0 

0 

0 

127 

21 . 8754 

200.5049 

36 

G08 

1 

10 

0 

0 

0 

101 

21.8662 

200 . 5413 

29 

G08 

1 

10 

0 

0 

0 

129 

21 .8698 

200.5786 

32 

GO  9 

1 

330 

0 

0 

0 

131 

21 . 9520 

200 . 6557 

10 

NA 

1 

315 

0 

0 

0 

where  NrSta  is  the  name,  DepNr  the  water  depth,  LatN  the  latitude  in 
degrees  North,  and  LonE  the  longitude  in  degrees  east  of  the  nearshore 
station.  OffSta  is  the  name  of  offshore  station  (e.g.,  WISWAVE  output  grid 
point).  The  variable  Opt  is  used  to  define  the  input  wave  parameters  for 
WAVETRAN.  Table  5  shows  the  options  for  the  parameter  Opt  and 
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Table  5.  Options  for  input  wave 
parameters  for  WAVETRAN. 


corresponding  wave  input  parameters  that  must  be 
specified.  The  subscripts  s,  p,  and  m  in  Table  5  refer 
to  the  significant,  peak,  and  mean  wave  parameters, 
respectively. 

The  orientation  of  offshore  bathymetric  contours  is 
specified  with  the  angle  AngBath  (slope  angle),  defined  as  the  angle  of  the 
up-slope  direction  of  the  bathymetry  measured  counterclockwise  from 
North.  This  is  illustrated  in  Figure  7. 


Opt 

Wave  Parameters 

1 

Hs,  Tp,  Dirp 

2 

Hs,  Tp,  Dirm 

3 

Hs,  Tm,  Dirp 

4 

Hs,  Tm,  Dirm 

Figure  7.  Definition  of  AngBath  parameter  for  offshore 
bathymetry  orientation. 


The  variable  KSHN  is  used  to  define  wave  directions  that  are  sheltered  by 
islands  or  protruding  land  masses,  where  the  last  digit  may  assume  values 
of  N  =  o,  1,  2.  Therefore,  KSHo  =  o  is  for  no  sheltering,  1  for  one-sided 
sheltering,  and  2  for  two-sided  sheltering.  KSHi  defines  the  sheltered 
angles  to  an  observer’s  left  when  facing  offshore.  KSH2  defines  the  shel¬ 
tered  angles  to  an  observer’s  right,  both  with  respect  to  the  shoreline 
orientation.  Table  6  shows  the  sheltering  values  of  KSHN  for  different 
sheltering  angles.  The  wave  angles  used  for  sheltering  are  specified  with 
respect  to  the  coastline.  Figure  8  shows  an  illustration  of  the  wave  angle 
convention  for  wave  sheltering. 
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Table  6.  Wave  sheltering  options. 


Sheltering  Angle  Specification 

Sheltering 

1-Sided 

2-Sided 

Sheltered 

Option 

Sheltering 

Sheltering 

Wave  Angle 

(KSH) 

Key  (KSH  1) 

Key  (KSH  1) 

(deg) 

0 

0 

0 

None 

1 

1 

0 

0-10 

1 

2 

0 

0-20 

1 

3 

0 

0-30 

1 

4 

0 

0-40 

1 

5 

0 

0-50 

1 

6 

0 

0-60 

1 

7 

0 

0-70 

1 

8 

0 

0-80 

1 

9 

0 

0-90 

1 

10 

0 

80-180 

1 

11 

0 

90-180 

1 

12 

0 

100-180 

1 

13 

0 

110-180 

1 

14 

0 

120-180 

1 

15 

0 

130-180 

1 

16 

0 

140-180 

1 

17 

0 

150-180 

1 

18 

0 

160-180 

1 

19 

0 

170-180 

2 

1  -  9 

10 

KSH1  value +  80-  180 

2 

1  -  9 

11 

KSH1  value +  90-  180 

2 

1  -  9 

12 

KSH1  value  +  100-  180 

2 

1  -  9 

13 

KSH1  value +  110 -180 

2 

1  -  9 

14 

KSH1  value  +120-180 

2 

1  -  9 

15 

KSH1  value  +  130-  180 

2 

1  -  9 

16 

KSH1  value  +  140-  180 

2 

1  -  9 

17 

KSH1  value  +  150-  180 

2 

1  -  9 

18 

KSH1  value  +160-  180 

2 

1  -  9 

19 

KSH1  value +  170 -180 

Figure  8.  Wave  angle  convention  for  wave  sheltering. 
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Transects  file 

The  wave  transformation  and  calculation  of  ponding  levels  over  reefs  and 
runup  require  the  specification  of  the  following  parameters:  shoreline 
orientation  ( AngSho ),  water  depth  over  the  reef  (DepRf),  width  of  the  reef 
(WthRf),  diffraction  coefficient  ( Diff ),  roughness  factor  (Rough)  for  runup 
calculations,  and  a  factor  that  accounts  for  the  berm  effects  (Berm).  The 
default  values  of  Diff,  Rough,  and  Berm  are  l.  The  convention  for  the 
shoreline  orientation  is  the  same  as  for  the  offshore  bathymetric  contours 
(Figure  7).  Guidance  on  estimates  of  the  roughness  factor  and  the  berm 
factors  is  available  in  the  Coastal  Engineering  Manual  (HQUSACE  2006) 
Sections  \T-5-2-a  and  5-3-a.  The  wave  conditions  for  nearshore  wave 
transformation  and  selected  transects  are  specified  with  the  variable 
StaNr.  An  example  of  the  Profile  Information  File  is  provided  next. 


Profile 

AngSho 

StaNr 

WthRf 

DepRf 

Slope 

Diff 

Rough 

Berm 

P030 

5.0 

130 

0 

0 

0 

1 

1 

1 

P003 

11.0 

103 

0 

0 

0 

1 

1 

1 

P028 

7 . 0 

128 

0 

0 

0 

1 

1 

1 

P027 

5.0 

127 

0 

0 

0 

1 

1 

1 

P001 

15.0 

101 

0 

0 

0 

1 

1 

1 

P02  9 

328.0 

129 

0 

0 

0 

1 

1 

1 

The  Profile  Information  File  should  be  located  in  the  input  folder  of 
TWAVE.  Similar  to  the  Station  Information  File,  the  Transect  File  is  also 
saved  in  the  Excel  worksheet  (labeled  with  the  variable  PROF)  for 
reference  and  visualization  purposes.  TWAVE  programs  cannot  read  the 
input  data  from  the  Excel  worksheets,  and  data  for  these  programs  must 
be  entered  in  the  corresponding  ASCII  files. 


Transects  bathymetry 

The  bathymetry  (topography)  for  each  transect  is  specified  in  the  Pi'ofile 
Information  File  using  a  prefix  and  file  extension  *.xy.  The  file  has  two 
data  columns,  and  the  user-specified  distances  are  in  the  first  column,  and 
elevations  are  in  the  second  column.  Distances  are  positive  in  the  offshore 
direction  as  measured  from  the  position  of  the  MSL  at  shoreline.  Eleva¬ 
tions  are  measured  with  respect  to  MSL  (below  the  MSL  values  are  nega¬ 
tive).  These  files  should  be  placed  in  the  folder  named  BoussiD. 


It  is  not  necessary  for  the  points  to  be  spaced  at  any  regular  interval 
because  numerical  wave  models  interpolate  transect  bathymetry  from 
these  files.  If  bathymetric  profile  data  for  a  new  transect  needs  to  be 


ERDC/CHL  TR-09-2 


42 


generated  from  raw  LIDAR  database,  it  may  be  necessary  to  apply  smooth¬ 
ing  to  the  extracted  profile  as  sharp  discontinuities  in  the  transect 
bathymetry  can  cause  instabilities  in  the  nearshore  wave  transformation 
models. 

BOUSS-1D  input  files 

As  mentioned  in  Chapter  4,  BOUSS-iD  is  run  for  each  transect  using  the 
values  of  nearshore  wave  height,  wave  period,  and  water  level  at  the  times 
of  maximum  runup  as  calculated  by  WAViD.  It  was  noted  in  Chapter  4 
that  it  is  unnecessary  to  run  BOUSS-iD  at  hourly  intervals.  Once  the  time 
periods  of  maximum  runup  are  determined  from  WAViD  results, 
BOUSS-iD  should  be  run  for  each  transect  at  these  times  with  the  wave 
height,  period,  and  water  levels  calculated  by  WAViD. 

BOUSS-iD  is  a  refined,  physics-based  advanced  wave  model  and  should  be 
run  when  the  maximum  runup  is  reached.  The  time  of  maximum  probable 
runup  is  determined  by  running  WAViD.  For  each  transect,  WAViD  saves 
the  wave  parameters  corresponding  to  maximum  probable  runup  in  a 
batch  file  named  as  xxx.bat.  This  batch  file  is  accessed  directly  by 
BOUSS-iD  through  a  setup  file.  The  BOUSS-iD  setup  file  is  the  only  input 
file  for  that  model.  It  contains  model  run  parameters,  post-processing 
tools  to  convert  model  binary  output  files  to  ASCII  files,  and  tools  for 
calculation  of  wave  runup  statistics.  Each  transect  has  its  own  setup  file 
named  as  Pppp_setup.gbat,  where  ppp  is  the  transect  ID  number.  An 
example  of  the  BOUSS-iD  script  file  is  shown  below. 
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(1)  #BOUSSlD_v3pO 

(2)  2  #  bathymetry  input  option  (1-2)  [1] 

(3)  %1  #  name  of  bathymetry  file  [.xy] 

(4)  %2  #  storm  surge/tidal  offset  (m)  [0.0] 

(5)  0.5  #  grid  spacing  for  numerical  computations  (m) 

(6)  -520.0  #  x  location  of  wave  generation  boundary 

(7)  2  #  type  of  wave  (1-2) 

(8)  #  time  series  synthesis  option  (1-2)  [1] 

(9)  #  incident  wave  spectra  option  (1-2)  [1] 

(10)  #  type  of  wave  spectrum  (1-5)  [1] 

(11)  #  JONSWAP  spectrum  option  (1-2)  [1] 

(12)  %3  #  peak  wave  period  (s) 

(13)  %4  #  significant  wave  height  (m) 

(14)  #  spectral  peakedness  parameter  gamma  [3.3] 

(15)  %3*0.6  #  minimum  wave  period 

(16)  #  maximum  wave  period  (s)  [25.0] 

(17)  #  Rescale  truncated  spectrum  (y/n) ?  [Yes]: 

(18)  300*%3  #  duration  of  synthesized  time  series 

(19)  #  time  step  (s) 

(20)  #  duration  of  numerical  simulation  (s) 

(21)  2.0  #  turbulent  length  scale  (m)  [0.09] 

(22)  20  #  Chezy  bottom  friction  factor  (10-1000) 

(23)  0.1  #  Smagorinsky  constant  [0.2] 

(24)  50  #  width  of  left  end  damping  layer  (m) 

(25)  #  damping  coefficient  for  left  end  damping  layer  (0-1)  [1] 

(26)  #  width  of  right  end  damping  layer  (m)  [0.0] 

(27)  0  #  damping  value  for  right  end  damping  layer  (0-1)  [1] 

(28)  #  number  of  instants  of  time  for  surface  elevation  output  [0] 

(29)  #  number  of  spatial  locations  for  time  series  output  [0] 

(30)  #  Create  an  animation  of  the  surface  elevation?  [No] : 

(31)  %1  #  prefix  for  output  files 

(32)  #END 

(33)  #EXPORT2 

(34)  %l_%5_xHsMWL.dat  #  name  of  ASCII  output  file  [.DAT] 

(35)  2  #  Implicit  Variable  Option  [1] 

(36)  2  #  number  of  GEDAP  input  files  (1-20)  [1] 

(37)  %l_hs  #  name  of  GEDAP  input  file  no.  1  [.001] 

(38)  %l_mwl  #  name  of  GEDAP  input  file  no.  1  [.001] 

(39)  #END 

(40)  #EXPORT2 

(41)  runup_ts.dat  #  name  of  ASCII  output  file  [.DAT] 

(42)  2  #  Implicit  Variable  Option  [1] 

The  first  column  (with  the  line  numbers  in  parentheses)  is  shown  here 
only  for  describing  the  lines,  but  this  column  does  not  exist  in  the  setup 
file  used  by  the  model.  Values  preceded  by  a  percent  sign  are  input 
parameters  that  are  passed  into  BOUSS-iD  model  from  the  batch  file  run 
boussid.bat.  This  batch  file  is  written  by  the  WAViD  model  and  the  user 
does  not  need  to  revise  the  file.  The  variables  passed  to  BOUSS-iD  from 
the  batch  file  are:  name  of  the  bathymetry  file  (line  3),  transect  label 
(JPppp,  line  31),  water  level  (line  4),  peak  wave  period  (line  12),  significant 
wave  height  (line  13),  and  storm  ID. 
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Lines  2  through  31  are  the  input  parameters  for  BOUSS-iD.  The  pound 
sign  (#)  in  front  of  a  program  name  signifies  executing  that  code.  The 
BOUSS-iD  model  input  parameters  specified  on  each  line  with  a  preceding 
pound  sign  (#)  are  user  comments  or  recommended  default  values.  A  brief 
description  of  input  parameters  is  provided  on  each  line  following  the 
pound  sign.  If  no  value  exists,  than  the  default  value  is  used  as  indicated 
by  the  values  in  the  square  brackets. 

The  values  of  significant  wave  height  and  mean  water  level  for  each 
transect  are  written  in  lines  33  through  38.  These  values  are  written  in  a 
file  named  Pppp_ssssss_xHsMWL.dat,  where  ppp  is  transect  ID  and 
ssssss  storm  ID.  Lines  39  through  43  are  saved  in  a  file  named 
runup_ts.dat  that  contains  wave  runup  time  series.  Line  44  executes  the 
Fortran  code  that  calculates  wave  runup  statistics.  For  additional  infor¬ 
mation,  see  the  Direct  Calculation  of  Wave  Runup  section  from  BOUSS-iD 
Results.  When  setting  up  a  new  transect,  it  is  recommended  to  run 
BOUSS-iD  interactively.  This  is  done  by  typing  the  executable  code  name 
in  a  DOS  window.  The  interactive  run  allows  the  user  to  physically  see  the 
values  of  various  parameters  as  calculated  by  BOUSS-iD  pre-processor. 
The  interactive  run  ensures  model  settings  are  correct,  including  the 
location  of  wavemaker,  width  of  the  damping  layers,  etc.  For  additional 
information  concerning  input  to  Boussinesq  model,  users  may  consult 
available  publications  (Demirbilek  et  al.  2005a  and  b;  Demirbilek  et  al. 
2007;  Demirbilek  and  Nwogu  2007;  Demirbilek  et  al.  2008).  Because 
BOUSS-iD  is  expected  to  be  used  for  many  transects  (tens  to  hundreds  of 
transects),  a  single  batch  file  called  run_bouss.bat  may  be  used  to  run 
BOUSS-iD  for  all  transects  at  once.  The  batch  file  allows  for  quick  turn¬ 
around  times  and  easier  editing  of  I/O  files. 

TWAVE  output  files 

TWAVE  results  can  be  viewed  and  plotted  using  Microsoft  Excel.  This 
eliminates  the  need  of  using  other  external  plotting  software.  Output  files 
are  imported  to  Excel  by  refreshing  the  external  data  read  by  worksheets 
and  selecting  the  appropriate  resource  files.  The  data  used  by  a  worksheet 
and  plots  are  also  automatically  updated.  Table  7  provides  a  description  of 
all  TWAVE  output  files,  including  the  file  names,  format,  location 
(directory),  and  associated  variables. 


Table  7.  TWAVE  output  file  names,  format,  location,  and  corresponding  TWAVE  variables. 


File(s) 

Description 

Default  Name  or  Format 

Location 

TWAVE 

Variable 

Last  Significant  Wave 
Height  File 

Contains  a  matrix  of  the  significant  wave  height  from  the  last 
time  step  of  WISWAVE  simulation 

Hs_last.dat 

.\offshore\ 

fhslast 

Maximum  Significant 
Wave  Height  File 

Contains  a  matrix  of  the  maximum  significant  wave  height 
during  the  WISWAVE  simulation 

Hs_max.dat 

.\offshore\ 

fhsmax 

Maximum  Wind  File 

Contains  a  matrix  of  the  maximum  wind  during  the  wind 
model  simulation 

Wmax.dat 

.\wind_fields\ 

fwmax 

Warm  Start  File 

Files  saved  and  read  by  WISWAVE  to  start  a  new  simulation 
at  an  intermediate  point  in  time 

ssssssmmddhh.mm 

.\offshore\ 

NA 

Offshore  Station  Time 
Series  Files 

Contains  a  time  series  of  wave  and  wind  conditions  at  the 
offshore  stations 

Gww_ssssss.off 

.\offshore\ 

NA 

Nearshore  Waves  File 

Contains  a  time  series  of  wave  conditions  at  individual 
nearshore  stations 

lxx_ssssss.PH3 

.\nearshore\ 

NA 

Nearshore  Wind  and 

Storm  Surge  Files 

Contains  a  time  series  of  wind,  atmospheric  pressure  and 
storm  surge  at  individual  nearshore  stations 

lxx_ssssss.uvs 

.\nearshore\ 

NA 

Nearshore  Astronomical 
Tides  File 

Contains  a  time  series  of  astronomical  tides  at  individual 
nearshore  stations 

lxx_ssssss.ast 

.\nearshore\ 

NA 

Nearshore  Station  Time 
Series 

Contains  a  time-series  of  the  waves,  wind,  storm  surge  and 
tide  at  nearshore  stations 

lxx_ssssss.wts 

.\nearshore\ 

NA 

Transect  Time  series 

Files 

Contains  a  time-series  of  the  wave  setup,  runup  and 
inundation  at  transects 

Pppp_ssssss.wsr 

.\transects\ 

NA 

Transect  Profile 

Contains  a  profile  of  the  wave  heights  and  water  levels  along 
each  transect  for  the  time  period  of  maximum  water  levels 

Pppp_ssssss_max.  dat 

.\transects\ 

NA 

WAV1D  Runup  Statistics 
File 

Contains  the  runup  statistics  calculated  by  WAV1D  at 
individual  transects 

Rstat_W  1  D_ssssss  .dat 

.\transects\ 

NA 

B0USS-1D  Transect  File 

Contains  the  cross-shore  variation  of  significant  wave  height 
and  setup  (mean  water  level) 

Pppp_ssssss_xHsmwl.dat 

.\BousslD\ 

NA 

B0USS-1D  Runup 
Statistics  File 

Contains  the  B0USS-1D  runup  statistics  for  all  profiles 

Rstat_B  1  D_ssssss  .dat 

.\BousslD\ 

NA 

Format  Key:  ssssss  -  storm  ID,  mm  -  month,  dd  -  day,  hh  -  hour,  ww  -  offshore  station  ID,  xx  -  nearshore  station  ID,  and  ppp  -  transect  ID. 

01 
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Deepwater  wave  and  wind  fields 

The  WISWAVE  model  calculates  values  of  the  maximum  significant  wave 
height  ( Hs,max )  and  significant  wave  height  at  the  last  time  step  ( Hs,iast )  at 
every  grid  cell.  These  results  are  saved  as  matrices  in  separate  ASCII  files 
in  the  folder  named  offshore.  The  default  names  of  these  output  files  are 
Hs_max.dat  and  Hs_last.dat,  respectively.  These  can  be  changed  in  the 
Control  File  using  the  variables  fhslast  and  fhsmax,  respectively.  The  grid 
cells  along  the  boundaries  and  cells  corresponding  to  land  are  assigned  a 
value  of  zero  in  the  files.  H s, max  and  Hs,iast  are  imported  to  the  worksheets 
labeled  HMAX  and  HLAST,  respectively.  The  PBL  wind  model  also  out¬ 
puts  the  maximum  wind  speed  at  every  grid  cell  during  the  simulation. 
This  file  is  saved  in  the  folder  named  wind_Jields,  and  the  default  file 
name  Wmax.dat  can  be  changed  in  the  Control  File  using  the  variable 
fwmax.  The  maximum  wind  speed  grid  can  be  imported  and  viewed  in 
Excel  using  the  worksheet  labeled  as  WMAX.  The  format  for  all  of  these 
files  is  shown  next,  where  X  represents  wave  heights  or  wind  speed  and 
the  number  in  parenthesis  are  the  values  of  cell  (I,J): 


X(l,  1) 

X (2, 1) 

X (3, 1) 

X  ( 4 , 1 )  . 

X  (n,  1 

X  ( 1 , 2  ) 

X (2, 2) 

X (3, 2) 

X  ( 4 , 2  )  . 

X  (n,  2 

X  ( 1 , 3 ) 

X (2, 3) 

X (3,  3) 

X  ( 4 , 3 )  . 

X  (n,  3 

X  ( 1 , 4  ) 

X (2, 4) 

X (3, 4) 

X  ( 4 , 4  )  . 

X  (n,  4 

X  ( 1 ,  m) 

X (2,m) 

X (3,m) 

X  ( 4 ,  m) 

X  (n,  m 

Offshore  stations  time  series 

This  file  contains  wave  information  at  offshore  stations  that  will  be  trans¬ 
formed  to  nearshore.  The  files  contain  values  of  the  wave  heights,  wave 
period,  and  wave  directions.  The  specific  format  and  variables  will  depend 
on  whether  the  WISWAVE  or  parametric  wave  model  is  run.  If  WISWAVE 
is  run,  the  post-processing  subroutine  spec__TW.exe  reads  the  WISWAVE 
output  file  and  writes  a  separate  output  file  for  each  offshore  station.  If  the 
parametric  wave  model  is  run,  these  files  are  used,  but  they  will  not  con¬ 
tain  wind  information.  These  files  are  saved  in  the  folder  offshore  and  the 
file  names  are  Gww_ssssss .  off,  where  ww  is  the  offshore  station  identifi¬ 
cation  number  and  ssssss  is  the  storm  identification  number.  An  example 
of  this  file  type  for  Hurricane  Iniki  is  provided  next. 
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199210  G26 
yyyymmddhh 

Hmo 

Tp 

Drp 

Tm 

Drm 

W 

Dir 

1992090919 

0.2 

4 

0 

4 

270 

7 

65 

1992090920 

0 . 4 

4 

54 

4 

281 

7 

65 

1992090921 

0.5 

4 

65 

4 

65 

7 

65 

1992090922 

0 . 6 

4 

65 

4 

65 

7 

65 

1992090923 

0.7 

4 

65 

4 

65 

8 

65 

1992091000 

CO 

o 

4 

65 

4 

65 

8 

65 

where  Hmo  =  zeroth  moment  wave  height  (m),  Tp  =  peak  period  (sec), 

Drp  =  peak  wave  direction  (degrees  clockwise  from  North),  Tm  =  mean 
wave  period,  Drm  =  mean  wave  direction  (degrees  clockwise  from  North), 
W  =  wind  speed  (m/sec),  and  Dir  =  wind  direction  (degrees  clockwise 
from  North).  The  first  line  indicates  the  storm  ID  followed  by  the  Offshore 
Station  ID.  The  second  line  of  the  file  is  a  header  line.  These  files  are 
imported  to  the  Excel  using  the  worksheet  labeled  OFF. 


Nearshore  waves  files 

These  files  contain  values  of  hourly  significant  wave  height,  wave  period, 
and  wave  direction  at  nearshore  stations  starting  at  the  first  hour  and 
ending  at  last  hour  of  the  Storm  Track  File,  respectively.  These  are  output 
by  WAVETRAN_TW.exe  model  that  performs  wave  transformation  from 
deep  water  to  shallow  water,  and  read  by  combine_TW.exe  model  that 
combines  winds,  waves,  tides,  and  storm  surge  information  into  a  single 
input  file  for  WAViD  to  use.  The  files  are  saved  in  the  folder  nearshore 
and  named  as  Ixx_ssssss.PH3,  where  xx  is  the  nearshore  station  ID,  and 
ssssss  is  the  storm  ID.  The  nearshore  wave  files  are  internal  files,  and 
users  need  not  be  concerned  with  these  files  unless  they  want  to  use  input 
wave  conditions  for  WAViD  from  a  different  source.  In  this  case,  users 
must  provide  a  file  for  each  nearshore  station  with  the  following  format: 


yyyymmddhh 

Hs 

T 

Theta 

1992091107 

0.6 

6.0 

131 . 4 

1992091108 

0.7 

6.0 

131 . 7 

1992091109 

00 

o 

6.0 

132.0 

1992091110 

0.9 

6.0 

132.3 

1992091111 

0.9 

7.0 

134 . 4 

1992091112 

1.0 

7.0 

135.0 

where  the  first  line  is  the  header  line,  followed  by  lines  containing  a 
10-digit  time  stamp,  the  significant  wave  height,  wave  period,  and  wave 
direction.  Here  and  in  the  subsequent  sections,  unless  otherwise  noted, 

T  =  wave  period  (sec),  Theta  =  wave  direction  (degrees  clockwise  from 
North).  Note  that  the  input  wave  height,  period,  and  direction  data  may  be 
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either  the  mean  or  peak  values  depending  on  the  opt  value  used  in  the 
deep  to  shallow  wave  transformation  (see  section  in  Chapter  4  called 
Deep-  to  Shallow-Water  Wave  Transformation).  Wave  parameters  are 
input  starting  1  hour  past  the  storm  track  time,  and  hourly  values  are 
required  through  the  end  of  the  storm  duration.  If  the  correct  time  stamp 
is  not  entered,  combine__TW.exe  will  give  an  error  message. 

Nearshore  wind  and  storm  surge  files 

These  files  contain  hourly  wind  speed,  direction,  atmospheric  pressure, 
and  storm  surge  at  individual  nearshore  stations  covering  the  same  time 
period  specified  in  the  Storm  Ti'ack  File.  These  files  are  located  in  the 
folder  nearshore  and  named  as  Ixx_ssssss.uvs,  where  xx  is  the  nearshore 
station  ID  and  ssssss  is  the  storm  ID.  These  files  are  internal  files  written 
by  the  wind  models  and  read  by  combine_TW.exe  model  to  write  the 
Nearshore  Stations  Time  Series  files  (see  the  section  Nearshore  Stations 
Time  Series  files  for  additional  information).  If  the  user  wants  to  enter 
different  wind  and  pressure  data  from  another  source  (model  or  database), 
this  is  done  by  writing  the  Nearshore  Wind  and  Storm  Surge  Files  in  the 
following  format: 


yyyymmddhh 

W 

Dir 

Pres 

Surge 

1992090918 

5.71 

22 . 8 

1009.8 

0.002 

1992090919 

5.85 

23.9 

1009.8 

0.002 

1992090920 

6.04 

25 . 1 

1009.8 

0.002 

1992090921 

6.28 

26.4 

1009.8 

0.002 

1992090922 

6.62 

27 . 8 

1009 . 7 

0.003 

where  the  first  line  is  the  header  line  that  provides  name  of  variables  for 
each  column,  Pres  =  atmospheric  pressure  (mbar),  Surge  =  storm  surge 
(m),  and  W  and  Dir  are  defined  as  in  previous  sections.  The  Nearshore 
Wind  and  Storm  Surge  files  are  internal  and  output  by  the  wind  models 
and  read  by  combine__TW.exe  model. 

Nearshore  astronomical  tide  files 

These  files  contain  time  series  of  the  astronomical  tides  at  each  nearshore 
station  at  hourly  intervals.  These  data  start  at  the  first  hour  and  end  at  the 
last  hour  specified  in  the  Storm  Ti'ack  File.  The  files  are  named  as 
Ixx_ssssss.ast  and  located  in  the  nearshore  folder.  The  program 
proc_tide.exe  converts  the  OTPS  output  file  to  a  Nearshore  Astronomical 
Tide  File.  The  program  proc_tide.exe  is  run  within  the  batch  script 
(run_otps.bat),  which  is  executed  from  the  Steering  Excel  worksheet 
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(STR).  The  batch  file  prepares  the  OTPS  input  files,  runs  the  OTPS  pro¬ 
gram,  and  runs  the  proc_tide.exe  code.  The  file  format  is  provided  below 
to  allow  users  to  implement  other  tidal  databases  that  can  be  written  in  the 
format  of  TWAVE  Astronomical  Tide  Files.  The  file  format  is: 


yyyymmddhh  Tide 

1992090918  -0.248 

1992090919  -0.235 

1992090920  -0.161 

1992090921  -0.040 


The  header  line  contains  the  time  stamp  and  the  variable  name  used  for 
astronomical  tides.  The  tides  should  be  in  m  and  referenced  to  the  local 
MSL. 

Nearshore  stations  time  series 

The  wave  conditions,  storm  surge  elevations,  and  tide  elevations  calcu¬ 
lated  from  different  components  of  TWAVE  are  combined  into  one  file  by 
the  program  combine_TW.exe,  and  these  files  are  saved  in  the  nearshore 
folder.  These  files  are  imported  to  the  Excel  worksheet  labeled  NR.  The 
format  name  for  these  files  is  Ixx_ssssss.wts,  where  xx  is  the  nearshore 
station  ID  and  ssssss  is  the  storm  ID.  The  top  section  of  the  Nearshore 
Station  time  series  file  is  shown  below  for  reference: 


199031  108 


yyyymmddhh 

H 

T 

Theta 

Tide 

W 

Dir 

Surge 

1990121918 

1.2 

11.0 

135.1 

-0 .468 

14 . 4 

181 . 6 

0 . 105 

1990121919 

1.2 

11 . 0 

136.4 

-0 . 391 

14 . 8 

180 . 4 

0 . 108 

1990121920 

1.2 

11 . 0 

136.4 

-0.264 

15.3 

179.1 

0 .111 

1990121921 

1.3 

11 . 0 

137 . 4 

-0 . 112 

15.8 

177 . 9 

0 .115 

1990121922 

1 . 4 

11.0 

137 . 4 

0 . 038 

16.3 

176.7 

0 .119 

The  first  line  contains  the  storm  ID  and  nearshore  station  ID  and  the 
second  line  is  a  header  line.  The  variable  names  are  defined  as  in  the 
previous  sections. 

Transect  time  series 

The  reef  wave  program__TW.exe  outputs  the  time  series  of  wave  setup, 
wave  runup,  and  inundation  at  each  transect  to  files  named  as 
~Pppp_ssssss.wsr,  where  ppp  is  the  transect  identification  number,  and 
ssssss  is  the  storm  ID.  The  files  are  saved  in  the  transects  folder  and  can 
be  imported  to  Excel  using  the  TRN  worksheet.  The  first  line  in  these  files 
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contains  the  storm  ID  and  the  transect  ID,  and  lines  2,  3,  5,  and  7  are  the 
header  lines. 

The  lines  8  and  beyond  contain  hourly  wave,  wind,  and  water  level  statis¬ 
tics  at  the  nearshore  station  and  over  the  reef-top,  and  runup  statistics. 

The  following  parameters  are  output  at  the  deepwater  end  of  transect: 
wave  height,  period,  direction,  wind  speed  and  direction,  tide,  and  storm 
surge.  The  following  parameters  are  output  over  the  reef-top:  wave  height 
Hr  (m),  setup  etar  (m),  and  water  depth  dr  (m).  The  2-percent  exceedance 
runup  {R.2.%)  is  in  meters  and  with  respect  to  the  SWL,  the  2-percent 
exceedance  water  level  {WL2%)  is  in  meters  and  with  respect  to  the  MSL. 
The  2-percent  exceedance  inundation  (72%)  is  in  meters  with  respect  to 
the  MSL  at  the  shoreline.  These  are  calculated  based  on  empirical  equa¬ 
tions  described  in  Chapter  4  using  the  wave  height  on  the  reef  top  (77r),  the 
deepwater  wave  length,  and  the  beach  slope. 

The  same  variables  as  in  the  hourly  time  series  are  output  for  the  highest 
WL%  and  R.2%,  and  these  are  shown  on  lines  4  and  6,  respectively.  An 
example  of  the  Transect  Time  Series  file  is  provided  next. 


199031  P001 


yyyymmddhh 

Hs 

Tp 

Theta 

w 

Dir 

Tide 

Surge 

etar 

dr 

Hr 

R2% 

WL2  % 

12% 

Maximum  Water  Level  Conditions 
1990122100  8.60  14.00  142.40  27.80 

306.20 

0.203 

0.232 

1.291 

1.527 

1.046 

1.701 

3.428 

74.08 

Maximum  Runup  Conditions 
1990122101  8.10  14.00  147.70 

26.10 

309.20 

0.255 

0.211 

1.190 

1 . 457 

0.996 

1 . 648 

3.304 

71.18 

Time-series 

1990122004 

3.80 

14.00 

116.20 

21.30 

171.30 

0.151 

0.161 

0.793 

0.906 

0.587 

1 . 173 

2.278 

45.41 

1990122005 

4 . 10 

14.00 

118.00 

22.40 

171.40 

0.113 

0.171 

0.843 

0.928 

0.607 

1.198 

2.325 

46.06 

1990122006 

4.50 

14.00 

119.40 

23.80 

171.70 

0.100 

0.184 

0.896 

0.981 

0.648 

1.249 

2.429 

47.52 

1990122007 

4.90 

14.00 

119.40 

25.20 

172.30 

0.116 

0.198 

0.940 

1.055 

0.701 

1.314 

2.568 

50.77 

1990122008 

5 . 10 

14 . 00 

118.00 

26.70 

173.00 

0.151 

0.214 

0.952 

1 .118 

0.742 

1.364 

2.681 

53.85 

1990122009 

5.30 

14 . 00 

116.20 

28.40 

174.40 

0.190 

0.232 

0.962 

1 .185 

0.787 

1 .416 

2.799 

57.23 

1990122010 

4.80 

14 . 00 

110.30 

30.30 

176.40 

0.212 

0.254 

0.879 

1 . 146 

0.752 

1.375 

2.719 

54 . 94 

1990122011 

4.80 

14.00 

108.80 

32.40 

179.30 

0.200 

0.279 

0.873 

1 . 153 

0.756 

1.379 

2.731 

55.27 

1990122012 

5.00 

14.00 

107.60 

34.70 

183.80 

0.144 

0.309 

0.903 

1 . 157 

0.763 

1.388 

2.743 

55.62 

1990122013 

5.30 

14.00 

107.60 

37.00 

190.30 

0.045 

0.340 

0.958 

1 . 144 

0.762 

1.387 

2.730 

55.24 

1990122014 

5.80 

14.00 

108.80 

39.70 

199.70 

-0.086 

0.377 

1.048 

1 .  140 

0.770 

1.396 

2.735 

55.38 

1990122015 

6.80 

14 . 00 

110.30 

42.00 

212.40 

-0.226 

0.408 

1.203 

1 .186 

0.814 

1 . 447 

2 . 832 

58 . 17 

1990122016 

7.80 

14 . 00 

113.20 

43.00 

228.70 

-0.349 

0.422 

1.360 

1.234 

0.856 

1 .495 

2 . 928 

60.93 

1990122017 

8.60 

14 . 00 

116.20 

42 . 10 

247.30 

-0.431 

0.411 

1 .488 

1.269 

0.886 

1.529 

2 . 997 

62.88 

1990122018 

9.10 

14.00 

119.40 

40.60 

262.90 

-0.455 

0.390 

1.565 

1.301 

0.910 

1.555 

3.055 

64.55 

1990122019 

9.00 

14 . 00 

123.30 

38.30 

275.20 

-0.414 

0.362 

1.547 

1.296 

0.906 

1.551 

3.046 

64.29 

1990122020 

9.00 

14.00 

130.60 

35.70 

285.90 

-0.316 

0.330 

1.511 

1.326 

0.925 

1.571 

3.096 

65.73 

1990122021 

8.90 

14.00 

136.50 

33.70 

293.10 

-0.180 

0.300 

1.443 

1.364 

0.947 

1.596 

3.159 

67.40 
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Transect  profile 

The  variations  of  wave  parameters  and  water  levels  along  each  transect  are 
output  for  the  hour  of  maximum  WL2%  to  files  named 
Pppp_ssssss_max.dat,  where  ppp  is  the  transect  ID.  These  files  are 
written  to  the  folder  named  transects  and  can  be  imported  and  viewed  in 
Excel  using  the  worksheet  named  PRO.  An  example  of  a  Transect  Profile 
File  is  shown  below. 


199031  P001 

yyyymmddhh  Hs  T  Theta  W  Dir  Tide+Surge 

1990122100  8.60  14.00  142.40  27.80  306.20  0.435 


xr 

hr 

dr 

Hr 

etar 

ar 

62  . 

99  -0.12 

1.697 

0.9414 

1.385 

-4.029 

R2  % 

Rmax 

WL2  % 

WLmax 

12% 

Imax 

1 . 07 

1.28 

2.89 

3.10 

108.39 

114.11 

X 

h 

slope 

d 

Hs 

eta 

a 

Qb 

L 

V 

-334 .01 

19.270 

0.318 

19.265 

8 

.  648 

-0.005 

-17 . 023 

0.091 

181 . 056 

0.248 

-333.51 

19 . 112 

0.318 

19.105 

8 

.  656 

-0.006 

-16.961 

0.096 

180.414 

0.256 

-333.01 

18 . 953 

0.318 

18 . 945 

8 

.  664 

-0.007 

-16.899 

0.102 

179.768 

0.265 

-332.51 

18.794 

0.318 

18.785 

8 

.  673 

-0.009 

-16.836 

0.108 

179.117 

0.274 

-332.01 

18 . 635 

0.318 

18  .  625 

8 

.  681 

-0.010 

-16.773 

0.114 

178 .462 

0.283 

-331.51 

18.476 

0.318 

18 .465 

8 

.  689 

-0.011 

-16.709 

0.120 

177 .802 

0.292 

The  storm  ID  number  and  transect  ID  are  in  the  first  line.  The  second  line 
is  the  header  line,  and  provides  names  of  the  nearshore  parameters  at  the 
hour  of  maximum  WL2%.  The  values  of  these  parameters  are  given  on  line 
three.  The  forth  line  is  the  header  line  for  the  reef-top  parameters  whose 
values  are  provided  on  line  five.  These  parameters  are:  xr  =  cross-shore 
location  where  the  reef-top  is  defined  (m),  hr  =  reef-top  water  depth 
(m,  MSL),  dr  =  reef-top  water  depth  with  respect  to  SWL,  Hr  =  wave  height 
over  the  reef-top  (m),  etar  =  wind  and  wave  setup  over  the  reef-top  (m), 
ar  =  reef-top  incident  wave  angle  (degrees  with  respect  to  shore  normal). 
The  sixth  line  is  the  header  line  for  the  runup  statistics,  and  values  of  these 
are  shown  on  line  seven.  The  eighth  line  is  the  header  line,  where  eta  = 
wave  +  wind  setup,  a  =  incident  wave  angle,  Qb  =  breaking  index  or 
fraction  of  broken  waves,  L  =  wave  length,  V  =  alongshore  current  speed. 

WAV1D  Runup  Statistics  File 

Wave  runup  statistics  are  calculated  in  WAViD  using  empirical  formulas 
described  in  the  section  “Calculation  of  Wave  Runup  Statistics,  Runup 
Time  Series.”  The  WAViD  Runup  Statistics  File  is  located  in  transect 
folder  and  may  be  viewed  in  Excel  using  the  worksheet  named  RWiD.  The 
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file  contains  the  transect  name  and  estimates  of  the  maximum  and 
2-percent  exceedance  runup,  water  levels,  and  inundation  limits. 


Trans 

Rmax 

R2% 

WLmax 

WL2  % 

Iraax 

12% 

P030 

4 . 81 

6.21 

6.18 

7.59 

36.05 

44.25 

P027 

5.96 

7.49 

6.70 

8.23 

61 . 64 

65.19 

P001 

3.32 

4.20 

4.49 

5.37 

43.01 

46.89 

B0USS-1D  transect  files 

The  BOUSS-iD  transect  files  contain  calculated  cross-shore  variation  of 
wave  heights  and  setup  for  each  transect.  The  file  names  are 
Pppp_ssssss_xHsmwl.dat,  where  ppp  is  transect  ID  and  ssssss  the  storm 
ID.  The  files  are  located  in  the  BoussiD  folder  and  may  be  imported  and 
viewed  in  the  Excel  worksheet  BiD.  These  ASCII  files  contain  three 
columns  of  data:  (l)  x-coordinate  (m),  (2)  zero-th  moment  wave  height 
(m),  and  (3)  total  mean  water  level  (MWL)  (m).  Note  the  wind  setup  is 
added  to  the  entire  domain  of  BOUSS-iD  and  MWL  has  the  wave  setup 
(e.g.,  MWL=SWL+wave  and  wind  setup). 


BOUSS-ID  Runup  Statistics  File 

Runup  statistics  are  calculated  from  BOUSS-iD  runup  time  series  files 
following  the  procedure  outlined  in  section  “Calculation  of  Wave  Runup 
Statistics,  Runup  Time  Series.”  The  BOUSS-iD  Runup  Statistics  File  is 
located  in  the  BoussiD  folder  and  may  be  imported  and  viewed  in  Excel 
using  the  worksheet  RBiD.  The  file  contains  the  transect  name,  maximum 
runup,  and  estimates  of  the  50-,  33-,  10-,  and  2-percent  exceedance  level 
runup.  An  example  BOUSS-iD  Runup  Statistics  File  is  shown  next. 
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6  Example  Applications 

Example  1:  Hurricane  Iniki  (Kauai) 

Hurricane  Iniki  is  the  most  recent  hurricane  that  has  made  landfall  on  the 
Hawaiian  Islands  and  produced  significant  flooding  and  wind  damage  on 
the  island  of  Kauai.  The  storm  was  upgraded  to  hurricane  strength  on 
9  September  1992  at  12:00  GMT  approximately  800  km  offshore.  It  made 
landfall  at  approximately  00:00  GMT  on  12  September  1992  as  a  cate¬ 
gory  4  storm  in  Saffir-Simpson  scale.  This  hurricane  produced  maximum 
winds  in  excess  of  120  knots.  Figure  9  shows  the  storm  track  for  Hurricane 
Iniki  and  locations  of  the  National  Data  Buoy  Center  (NDBC)  buoys  in  the 
area. 

In  the  Northern  Hemisphere,  the  strongest  winds  occur  in  the  forward 
right  quadrate  of  the  hurricane.  Thus,  the  southeast  portion  of  the  island 
received  the  strongest  winds  and  highest  waves.  Figure  10  shows  the  spa¬ 
tial  field  of  maximum  winds  at  each  grid  cell  as  computed  with  the  PBL 
model  from  hourly  values.  The  curved  contours  are  due  to  interpolation 
between  the  hourly  intervals. 

Wind  and  wave  data  were  collected  by  four  buoys  (Figure  9).  However, 
buoys  51001  and  51004  are  too  far  from  the  storm  track  and  could  not  be 
used  in  validation;  only  buoys  51002  and  51003  were  used  in  the  modeling 
validation.  The  buoys  recorded  hourly  wind  speed  and  direction  for 
8.5  min.  Wind  speeds  were  converted  from  the  buoy  anemometer  height  of 
5  to  10  m  using  the  method  of  Liu  et  al.  (1979).  Significant  wave  height  and 
mean  wave  period  were  also  recorded  at  hourly  intervals  as  20-min 
records. 
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Figure  9.  Storm  track  for  Hurricane  Iniki. 
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Figure  10.  Maximum  wind  speeds  for  Hurricane  Iniki  calculated  using  PBL  model. 


Figures  n  and  12  show  a  comparison  of  measured  and  computed  wind 
speeds  at  buoys  51002  and  51003  using  the  PBL,  Holland,  and  Modified 
Rankine  Vortex  models.  The  Holland  model  was  run  with  a  scaling 
parameter  of  B  =  1.3,  and  the  shape  parameter  for  the  Modified  Rankine 
Vortex  model  was  set  to  X  =  0.5. 

The  storm  passed  within  approximately  250  km  to  the  south  of 
Buoy  51002  and  460  km  to  the  east  of  Buoy  51003.  All  three  models 
performed  reasonably  well  compared  to  measured  wind  speeds  reported 
by  both  buoys.  The  Modified  Rankine  Vortex  and  Holland  Models 
produced  similar  wind  speeds  at  Buoy  51002  and  the  Holland  and  PBL 
models  produced  similar  wind  speeds  at  Buoy  51003.  The  wind  direction 
for  Modified  Rankine  Vortex  and  Holland  models  are  identical  because 
they  were  calculated  based  on  the  position  of  storm  relative  to  the  buoys 
(see  Chapter  2  for  details). 
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Figure  11.  Comparison  of  (a)  wind  speed  and  (b)  direction  at  Buoy  51002  for  various  wind  models. 
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Wind  Speed  at  Buoy  51003 
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Figure  12.  Comparison  of  (a)  wind  speed  and  (b)  direction  at  Buoy  51003  for  various  wind  models. 


WISWAVE  was  run  for  the  Hawaiian  Islands  with  the  PBL  model  winds. 
The  maximum  significant  wave  height  for  each  model  grid  point  is  shown 
in  Figure  13.  Calculated  significant  wave  heights  for  Hurricane  Iniki  from 
WISWAVE  were  greater  than  15  m  about  280  km  south  of  Kauai.  Near  the 
south  coast  of  Kauai,  significant  wave  heights  reached  13  m,  and  9  to  10  m 
on  the  east  and  west  coasts. 
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Maximum  Significant  Wave  Height,  m 
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Figure  13.  Maximum  significant  wave  heights  during  Hurricane  Iniki  calculated  using 

WISWAVE  and  PBL  winds. 


The  comparison  of  measured  and  computed  wave  heights  and  periods  are 
shown  in  Figures  14  and  15.  Significant  wave  height  is  defined  as  the  zero- 
moment  wave  height  here.  The  comparison  at  Buoy  51002,  which  was  to 
the  right  of  the  forward  direction  of  storm,  is  good.  The  computed  mean 
wave  periods  are  also  within  2  sec  of  the  measured  values.  The  computed 
wave  heights  compare  well  with  measurements.  The  maximum  difference 
between  calculated  and  measured  wave  heights  is  2m,  occurring  between 
0:00  and  12:  00  GMT  on  11  September  1992. 
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Figure  14.  Comparison  of  significant  wave  heights  and  mean  wave  periods  at  Buoy  51002 

for  Hurricane  Iniki. 
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Figure  15.  Comparison  of  significant  wave  heights  and  mean  wave  periods  at  Buoy  51003 

for  Hurricane  Iniki. 

Without  calibrating  models  to  data  from  this  storm,  a  much  greater  differ¬ 
ence  occurs  between  the  computed  and  measured  wave  heights  and  wave 
periods  to  the  right  of  the  storm.  The  computed  wave  heights  were  under¬ 
estimated  by  2m  before  io:oo  GMT  on  n  September,  and  over  estimated 
by  3m  when  the  storm  passed  close  to  the  buoy  at  around  07:00  GMT  on 
11  September  1992.  The  mean  wave  periods  have  a  maximum  difference  of 
approximately  4  sec.  The  under-  and  over-predicted  wave  height  differ¬ 
ences  are  significant.  There  was  also  a  phase  lag  that  caused  an  associated 
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lag  in  the  surge  peak.  Additional  validation  is  necessary  and  planned  to 
investigate  the  causes  of  these  large  differences.  However,  for  the  intended 
purposes  of  TWAVE  modeling  system,  these  non-calibrated  TWAVE 
modeling  results  are  comparable  to  those  reported  by  previous  modeling 
studies  (Phadke  et  al.  2003;  Cheung  et  al.  2003). 

The  wave  and  wind  parameters  are  plotted  in  Figure  16  at  the  offshore 
sta  G08,  located  just  south  of  the  Poipu  Beach.  Poipu  Beach  is  located  on 
the  southeast  coast  of  Kauai.  Poipu  Beach  received  some  of  the  largest 
waves  during  Hurricane  Iniki  because  of  its  location  with  respect  to  the 
storm  track.  The  maximum  wind  and  wave  heights  occurred  during  the 
landfall  at  00:00  GMT  on  12  September  1992. 

As  the  storm  makes  landfall,  the  wave  and  wind  directions  change  sharply. 
The  mean  wave  direction  is  a  better  parameter  than  the  peak  wave  direc¬ 
tion  for  nearshore  wave  transformation  as  it  appears  to  be  more  stable  and 
follows  closely  the  wind  direction.  Figure  17a  shows  the  nearshore  water 
levels  at  the  sta  I01  which  is  located  south  of  Poipu  Beach.  The  astronomi¬ 
cal  tide  was  calculated  using  the  Oregon  State  Tidal  Prediction  Software 
and  the  Hawaii  regional  tidal  database.  The  storm  surge  was  calculated 
using  only  the  atmospheric  pressure  (barometric  tide).  The  storm  landfall 
occurred  shortly  before  the  high  tide  stage.  The  offshore  wave  conditions 
were  transformed  to  the  nearshore  using  WAVETRAN  (Gravens  et  al. 
1991).  The  nearshore  wave  height  and  wave  period  are  shown  in 
Figure  17b. 

The  maximum  wave  heights  occurred  on  00:00  GMT  on  12  September 
1992  at  the  offshore  station.  Because  of  the  influence  of  the  wave  direction, 
the  maximum  wave  height  at  the  nearshore  station  occurred  1  hr  later. 
Most  of  the  wave  energy  is  concentrated  to  within  a  16-hr  period  centered 
on  the  landfall  time  of  storm.  Port  Allen  is  located  in  Hanapepe  Bay  on  the 
south  central  coast  of  Kauai.  Water  levels  were  recorded  in  Port  Allen  at 
NOAA  tide  gauge  sta  1611347.  The  measured  and  calculated  water  levels  at 
Port  Allen  are  shown  in  Figure  18.  The  Total  water  level  is  the  sum  of  the 
astronomical  tide,  storm  surge,  and  setup  (wind  plus  wave  setup).  The 
computed  and  measured  total  water  levels  show  a  good  agreement.  The 
time  shift  in  the  peak  water  levels  might  be  due  to  the  simplified  storm 
surge  calculation  used  in  this  example. 
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Figure  16.  Offshore  sta  G08  during  Hurricane  Iniki.  (a)  Zeroth  moment  wave  heights 
(Hmo),  peak  wave  periods  (TP),  and  mean  wave  period  (Tm),  (b)  wind  speed, 
and  (c)  peak  wave  period  (Drp),  mean  wave  period  (Dm),  and  wind  direction  (Dir). 
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Figure  17.  Tide,  surge  and  combined  surge,  and  tide  water  elevations:  (a)  wave  heights 
and  periods,  (b)  at  nearshore  sta  101. 
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Figure  18.  Comparison  of  measured  and  computed  water  levels  and  Port  Allen,  Kauai 

during  Hurricane  Iniki. 
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Because  Poipu  Beach  is  located  on  the  southeast  coast  of  Kauai,  it  received 
some  of  the  largest  waves  during  Hurricane  Iniki.  A  time  series  of  water 
levels,  wave  parameters,  and  2-percent  exceedance  water  levels  (fVL2%)  at 
Poipu  Beach  are  shown  in  Figure  19.  Here  variables  dr  and  Hr  are  the 
water  depth  and  wave  height  over  the  reef  top.  Figure  19  shows  the  maxi¬ 
mum  coastal  inundation  occurred  during  storm  landfall  at  02:00  GMT  on 
12  September  1992. 
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Figure  19.  Time  series  statistics  at  Poipu  Beach  during  Hurricane  Iniki. 

The  WAViD  model  offers  an  efficient  method  for  obtaining  the  time  period 
of  maximum  coastal  inundation.  This  significantly  reduces  the  number  of 
simulations  required  using  more  computationally  intensive  lD  or  2D 
models.  Although  the  maximum  storm  surge  and  wave  setup  occurred 
near  high  tide  for  Hurricane  Iniki,  this  is  not  always  the  case  for  storm 
events.  The  wave  runup  is  a  nonlinear  function  of  water  levels,  wave 
heights,  and  wave  periods.  By  using  simple  empirical  formulas,  WAViD 
can  predict  runup  and  inundation  for  the  entire  storm  duration  and  a 
reasonably  accurate  estimate  of  time  for  the  highest  coastal  flooding  to 


occur. 
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The  cross-shore  variation  of  wave  height  (Hs),  wave  breaking  (Qi,),  and 
setup  (eta)  at  Poipu  Beach  as  calculated  from  WAViD  are  shown  in 
Figure  20a.  The  bathymetric  profile  is  shown  in  Figure  20b.  Here  Qb  is  to 
the  fraction  of  breaking  waves  whose  value  varies  between  o  and  1.  Near¬ 
shore  waves  propagating  over  the  reef  shoal  and  increase  in  height,  and 
eventually  the  breaker  zone  starts  at  about  400  m  from  the  shoreline.  The 
wave  height  and  setup  over  the  reef  as  predicted  from  WAViD  are  approx¬ 
imately  1.5  and  0.6  m,  respectively.  The  maximum  water  level  predicted  by 
WAViD  using  the  empirical  formula  of  Mase  (1989)  is  5.5  m,  which  com¬ 
pares  well  to  the  measured  maximum  runup  from  the  debris  line  of  4-6  m 
reported  by  Fletcher  et  al.  (1995). 


Distance,  m 


Figure  20.  Wave  transformation  over  Poipu  reef  for  Hurricane  Iniki  calculated  from  WAVID 

on  12  September  1992,  00:00  GMT. 
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The  nearshore  transformation  was  also  simulated  using  BOUSS-iD  for  the 
(wave  and  water  level  conditions  at  the  time  of  maximum  coastal  flooding 
as  predicted  by  WAViD.  The  calculated  cross-shore  variation  of  wave 
height  (Hs)  and  setup  (eta)  are  shown  in  Figure  21.  The  results  from 
WAViD  and  BOUSS-iD  models  are  relatively  similar.  One  difference  is 
that  the  BOUSS-iD  model  does  not  predict  an  increase  in  wave  height 
before  breaking.  This  may  be  due  to  the  values  of  model  parameters  used 
in  these  simulations.  The  predicted  wave  height  and  setup  over  the  reef 
from  BOUSS-iD  are  1.5  and  0.4  m,  respectively,  and  similar  to  those 
estimated  by  WAViD  model. 
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Figure  21.  Wave  transformation  over  Poipu  reef  for  Hurricane  Iniki  calculated  from  BOUSS-ID 

on  12  September  1992,  00:00  GMT. 
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Example  2:  Typhoon  Russ  (Guam) 

Typhoon  Russ  passed  approximately  70  km  south  coast  of  the  island  of 
Guam  on  20  December  1990,  moving  from  east  to  west  and  with  maxi¬ 
mum  l-min  sustained  wind  speeds  of  125  knots.  This  powerful  storm 
produced  sustained  wind  speeds  in  excess  of  100  knots  on  Guam  during  its 
passage  along  a  track  of  about  35  miles  south  of  the  U.S.  Territory  of  Guam 
(Figure  22).  Modeling  results  are  shown  in  Figure  23  at  sta  I01  for  the 
section  of  Guam  coast  that  lies  between  Inarajan  Bay  and  Agfayan  Bay  in 
the  southeastern  part  of  the  island.  This  section  of  coast  is  directly  exposed 
to  typhoons  approaching  the  island  from  the  east.  The  coastal  roads  and 
properties  in  this  area  experience  damage  from  typhoon  inundation. 


Storm  Track  -  Typhoon  Russ 
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Figure  22.  Storm  track  for  Typhoon  Russ  (red  line  indicates  coastline  for  Guam,  green  line 

indicates  coastline  for  Rota). 
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Nearshore  Stations 


Longitude,  °E 

Figure  23.  Island  of  Guam  showing  location  of  nearshore  stations 

Figure  24  shows  the  contour  map  of  the  maximum  significant  wave 
heights  at  each  WISWAVE  grid  cell.  This  figure  is  plotted  using  the  Excel 
worksheet  HMAX.  The  maximum  calculated  significant  wave  heights 
reached  18  m  in  the  southeast  of  Guam,  and  7.5-10  m  close  to  the  shore. 

Figure  25  shows  the  time  series  of  wave  heights,  periods  and  directions  at 
the  offshore  sta  G26.  This  station  is  located  on  the  southeast  side  of  the 
island  of  Guam  in  the  Inarajan  Bay.  The  maximum  calculated  offshore 
significant  wave  heights  at  G26  are  15  m,  and  the  peak  wave  period  and 
direction  are  12  sec  and  160  deg,  respectively. 
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Figure  24.  Maximum  significant  wave  heights  during  Typhoon  Russ  near  the  island  of  Guam. 


Hourly  wave  and  water  level  conditions  are  plotted  in  Figure  26  at  near¬ 
shore  sta  lot.  The  wave  conditions  were  calculated  using  a  spectral  wave 
transformation  from  the  offshore  sta  G26  to  the  nearshore  sta  lot  by 
assuming  straight  and  parallel  bathymetric  contours  (see  Deep-  to 
Shallow-Water  Transformation  section  in  Chapter  4).  With  the  type  of 
offshore  to  nearshore  transformation  performed,  the  significant  wave 
heights  were  reduced  from  15  to  9  m.  The  peak  storm  surge  coincided  with 
the  low  water  astronomical  tides  at  16:00  GMT  on  20  December  1990. 
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Figure  25.  Wave  conditions  at  Offshore  sta  G26:  (a)  Zeroth  moment  wave  heights  (Hmo),  mean 
wave  period  (Tm),  and  peak  wave  period  (TP),  (b)  peak  wave  direction  ( Drp )  and  mean  wave 

direction  (D  rm). 
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Figure  26.  Nearshore  sta  101  near  Inarajan  Guam  during  Typhoon  Russ:  (a)  wave  heights  and 
periods,  (b)  wave  direction,  and  (c)  tide,  surge  and  combined  surge  and  tide  water  elevations. 
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For  the  transect  Pool,  the  hourly  time  series  of  tides,  storm  surge,  reef  top 
water  depths,  significant  wave  height,  and  2-percent  exceedance  water 
levels  as  calculated  from  WAViD  are  shown  in  Figure  27.  This  figure 
illustrates  that  the  maximum  coastal  inundation  could  occur  at  times 
different  than  the  time  of  peak  storm  surge.  For  Typhoon  Russ,  the 
computed  storm  surge  peaked  at  16:00  GMT  on  20  December,  whereas 
the  calculated  maximum  coastal  inundation  occurred  6  hr  later. 


0:00  0:00  0:00  0:00  0:00  0:00  0:00  0:00 

Date,  GMT 


Figure  27.  Tide,  surge,  wave  setup  and  runup  at  transect  P001  near  Inarajan  Guam  during 

Typhoon  Russ. 


The  cross-shore  variation  of  the  significant  wave  height  and  mean  water 
level  for  transect  Pool  is  shown  in  Figure  28.  The  significant  wave  height 
and  mean  water  level  over  the  reef  are  approximately  1.0  and  1.3  m, 
respectively.  The  strongest  wave  breaking,  indicated  by  the  fraction  of 
breaking  waves  Qb,  occurs  at  the  reef  edge  (125  m  from  the  shoreline),  and 
at  the  beach  face,  but  wave  breaking  starts  as  far  as  350  m  from  the  shore 
at  approximately  15-m  water  depth. 
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Figure  28.  Transect  POOl  during  period  of  maximum  water  level  reached  by  runup  near 
Inarajan  Guam  during  Typhoon  Russ  on  12  December  1990,  00:00  GMT:  (a)  significant  wave 
height  Hs  (m),  setup  eta  (m)  and  fraction  of  breaking  waves  Qb  and  (b)  still  water  depths. 
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7  Summary  and  Conclusions 

A  user-friendly  modeling  package  called  TWAVE  is  described  in  this 
report.  This  modeling  system  is  developed  for  tropical  cyclones,  and  pro¬ 
vides  engineers  and  decision-makers  estimates  of  associated  winds,  waves, 
storm  surges,  and  resulting  coastal  inundation  levels.  TWAVE  consists  of  a 
family  of  computational  models,  and  the  system  is  implemented  in  the 
Microsoft  Excel  for  graphical  interface.  This  report  serves  as  the  User's 
Guide  and  presents  description  of  the  TWAVE  methodology,  input  files, 
output  files,  and  two  example  applications. 

The  TWAVE  system  provides  users  with  models  of  varying  complexity  and 
accuracy  for  calculating  winds,  waves,  and  coastal  inundation.  Different 
levels  of  complexity  are  provided  for  different  types  of  applications  such  as 
modeling  past  storm  events  (hindcast),  forecast  of  future  storms,  and 
hypothetical  storms.  Less  resource-demanding  models  that  run  fast  with  a 
low-order  accuracy  may  be  used  to  perform  quick  simulations  for  decision¬ 
making  and  determine  if  further  modeling  using  models  of  higher  accu¬ 
racy  is  warranted.  For  example,  TWAVE  parametric  models  may  be  used 
to  calculate  rough  estimates  of  coastal  inundation  at  hourly  intervals,  and 
the  approximate  time  of  the  most  severe  inundation  to  occur.  The  user 
may  then  choose  to  run  other  TWAVE  models  with  higher  accuracy 
(physics-based  models  with  fewer  free  parameters  that  are  resource¬ 
demanding)  to  pinpoint  the  locations  and  times  of  the  flooding  and  sever¬ 
ity  of  coastal  inundation.  This  approach  will  allow  users  to  reduce  the 
number  of  numerical  simulations  required  with  comparatively  more 
resource-demanding  numerical  models. 

Two  example  applications  presented  in  this  report  are  intended  to  illus¬ 
trate  the  usage  of  TWAVE  in  practical  applications.  Computed  values  of 
wave  heights,  water  levels,  wave  setup,  and  wave  runup  from  TWAVE  are 
used  to  calculate  the  extent  of  island  flooding,  and  compared  to  the  post¬ 
storm  damage  survey  records.  The  estimates  are  reasonable  in  spite  of 
several  uncertainties  involved,  and  demonstrate  the  usefulness  of  TWAVE 
for  quantitative  flooding  of  coastal  areas.  Guidance  is  provided  in  the 
example  applications  for  selecting  critical  modeling  parameters  including 
for  determining  the  breaking  wave  height  Hb  and  stable  wave  height 
parameter  used  in  WAViD  model.  In  the  next  update  of  TWAVE  model, 
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field  and  laboratory  data  will  be  used  to  calibrate  models  and  parameters 
for  breaker  types,  reef  characteristics  (reef  surface  irregularities  and 
roughness),  and  wind  conditions.  The  integrated  set  of  offshore  and  near¬ 
shore  models  for  winds,  tides,  waves,  and  flooding  as  implemented  in  the 
TWAVE  modeling  system  constitute  a  useful  predictive  engineering  tool 
for  planning  and  emergency  response  management  of  the  tropical  storms 
that  affect  islands  in  both  the  Pacific  Ocean  and  Caribbean  Sea.  TWAVE 
can  also  be  used  in  other  areas  of  the  U.S.  coasts  and  engineering  studies 
concerned  with  the  design  of  shore  and  wetland  protection  systems  such 
as  the  barrier  islands,  levees,  and  coastal  navigation  and  harbor  structures 

In  closing,  although  TWAVE  is  currently  configured  for  the  U.S.  Territory 
of  Guam  and  Kauai  Island,  HI,  it  can  be  easily  adapted  to  other  islands  or 
coastal  applications  for  prediction  and  risk  assessment  of  hurricane  and 
typhoon-induced  coastal  flooding.  Presently,  TWAVE  has  been  validated 
against  measured  winds,  waves,  and  coastal  inundation  records  obtained 
for  Hurricane  Iniki  that  made  landfall  on  the  island  of  Kauai,  HI,  in  the 
morning  of  12  September  1992.  Even  though  there  were  no  data  available 
for  Typhoon  Russ  that  impacted  the  U.S.  Territory  of  Guam  Island,  this 
storm  event  is  also  simulated  to  demonstrate  applicability  and  capabilities 
of  TWAVE  to  other  islands  and  storm  types  (typhoons).  For  Hurricane 
Iniki,  the  mean  wave  periods  as  compared  to  the  buoys  data  were  off 
approximately  4  sec,  wave  heights  differences  were  under-predicted  about 
2m  and  over-predicted  up  to  3m,  and  TWAVE  results  lagged  the  surge 
peak.  The  observed  lag  in  the  phase  gave  rise  to  an  associated  lag  in  the 
peak  surge  estimate.  Additional  validation  is  necessary  and  planned  to 
investigate  potential  causes  of  these  differences.  Overall,  the  non- 
calibrated  modeling  results  from  TWAVE  are  comparable  to  those 
reported  by  previous  modeling  studies. 

This  is  a  manual  for  the  TWAVE  model,  not  a  full  assessment  of  the 
model’s  abilities,  which  is  a  task  that  is  just  beginning,  now  that  the  model 
is  operational.  The  TWAVE  modeling  system  is  expected  to  evolve  and  be 
improved  as  it  is  used  in  engineering  projects.  TWAVE  will  be  calibrated 
with  additional  field  observations  in  progress.  Other  numerical  models 
and  different  engineering  methods  may  be  incorporated  into  TWAVE  as 
these  validations  warrant.  The  feedback  and  suggestions  from  the  users 
community  will  be  used  in  the  future  revisions  of  TWAVE.  Please  direct 
inquiries  and  suggestions  to  the  attention  of  authors. 
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